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. We study massive one-loop integrals by analytically continuing the Feynman inte- 

gral to negative dimensions as advocated by Halliday and Ricotta and developed by 
Qh| Suzuki and Schmidt. We consider n-point one-loop integrals with arbitrary powers of 

Q_i| propagators in general dimension D. For integrals with m mass scales and q external 

momentum scales, we construct a template solution valid for all n which allows us to 
obtain a representation of the graph in terms of a finite sum of generalised hypergeo- 
metric functions with m+q — 1 variables. All solutions for all possible kinematic regions 
are given simultaneously, allowing the investigation of different ranges of variation of 
mass and momentum scales. 

As a first step, we develop the general framework and apply it to massive bubble and 
vertex integrals. Of course many of these integrals are well known and we show that 
the known results are recovered. To give a concrete new result, we present expressions 
for the general vertex integral with one off-shell leg and two internal masses in terms of 
hypergeometric functions of two variables that converge in the appropriate kinematic 
regions. The kinematic singularity structure of this graph is sufficiently complex to 
give insight into how the negative-dimension method operates and gives some hope 
that more complicated graphs can also be evaluated. 
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1 Introduction 



Loop integrals play an important role in making precise perturbative predictions in quan- 
tum field theory, in general, and in the Standard Model of particle physics, in particular. 
As such, a large effort has been expended in developing methods for evaluating them. The 
problem is complicated by the appearance of ultraviolet (UV) and infrared (IR) singulari- 
ties, and it has become customary to use dimensional regularisation ^ ||, || to extend the 
dimensionality of the loop integral away from 4-dimensions to D = 4 — 2e, to regulate the 
infrared and ultraviolet singularities. 

With the increasing of the number of legs, of the number of mass scales or of the number 
of loops, the integrals can be made almost arbitrarily complex and difficult to solve analyt- 



ically. Different methods |4fl-[ll] have been developed to solve the Feynman integrals. We 



mention here only two of them: the integration by parts |7]], which works well for some two- 



loop vertex diagrams |L2| reducing them to simpler known graphs with different powers of 
the propagators (however this breaks down for more complicated graphs such as the double 
box, where irreducible numerators factors are present) and the Mellin-Barnes integral repre- 



sentation (see for example [0), which was successfully used by Smirnov ||T4| to calculate the 
two-loop box integral. In this approach, the integral is usually written as multiple contour 
integrals of T functions and powers of ratios of the mass scales in the problem. By closing 
the contour, we obtain an infinite series of residues at the singular points of the T functions. 
These series can be identified as generalised hypergeometric functions, whose convergence 
properties reflect the threshold-singularity structure of the integral. 

There are several advantages in using hypergeometric functions to represent the integral. 
First, these hypergeometric functions often have integral representations themselves, in which 
an expansion in e can be made, yielding expressions in logarithms, dilogarithms etc.. It seems 
that, where direct evaluation of the hypergeometric function in terms of known functions is 
possible, very compact results are obtained [ITH [RSI]. Second, because the series is convergent 



and well behaved in a particular region of phase space, it can be numerically evaluated JT7 . 
In fact, each hypergeometric representation immediately allows an asymptotic expansion 
of the integral in terms of ratios of momentum and mass scales. Third, through analytic 
continuation formulae, the hypergeometric functions valid in one kinematic domain can be 
re-expressed in a different kinematic region. 

Not all work has concentrated around D = 4. In fact, a close connection between tensor 
loop integrals - those with additional powers of the loop momentum in the numerator - and 



higher-dimension scalar integrals (D = 6 — 2e, for example) is well established [IT], [L8[ [19]. ^0 
Furthermore, in 1987, Halliday and Ricotta |2TJ, [22J suggested that it would be useful to 
calculate the loop integral considering D as a negative number. Because loop integrals are 
analytic in the number of dimensions D (and also in the powers of the propagators) they 
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proposed to calculate the integral in negative dimensions and return to positive dimensions, 
and specifically D = 4 — 2e, after the integrations have been performed. As we will discuss 
more fully later on, integration over the loop momentum and/or the parameters introduced 
to do the loop integration is replaced with infinite series, which again can be identified as 
generalised hypergeometric functions. Recently this idea has been picked up again by Suzuki 
and Schmidt who have evaluated a number of two- loop integrals |H§ , three-loop integrals pi 



one-loop tensor integrals [25] as well as the one- loop massive box integral for the scattering 



of light by light [26|. In this latter well as reproducing the known hypergeometric- 



series representations of Ref. [fL5 |, valid in particular kinematic regions, Suzuki and Schmidt 



simultaneously found hypergeometric solutions valid in other kinematic domains. Of course, 
all of these solutions are related by analytic continuation. However, it is easy to envisage 
integrals that yield hypergeometric functions where the analytic continuation formulae are 
not known a priori. In these cases, having series expansions directly available in all kinematic 
regions is useful. 

In this paper we wish to explore the negative-dimension approach (NDIM) further. In 
particular we focus on one-loop integrals with general powers of the propagators and arbi- 
trary dimension D. There are several reasons for doing this. First, it allows connection with 
the general tensor- reduction program based on integration by parts of Refs. fl9| , 20|| . Here 



the tensor integrals are linear combinations of scalar integrals with either higher dimension 
or propagators raised to higher powers. Second, we can imagine inserting the one-loop re- 
sults into a two-loop integral by closing up external legs. This is trivial for most bubble 
integrals, but more complicated for vertex and box graphs. Broadhurst [27J has shown that 



this is possible for the non-trivial two-loop self-energy graph. Third, it actually simplifies 
the calculation. As we will show, by keeping the parameters general, it is easier to identify 
the regions of convergence of the hypergeometric series and therefore which hypergeometric 
functions to group together. For specific values of the parameters, the hypergeometric func- 
tions often collapse to simpler functions. As a first step, we develop the general framework 
and apply it to massive bubble and vertex integrals. Of course many of these results are well 
known. However, they serve to iron out some of the subtleties of the NDIM approach. To 
give a concrete new result, we present expressions for the vertex integral with one off- shell 
leg and two internal masses in arbitrary dimension and for general powers of the propagators 
in terms of hypergeometric functions of two variables. The kinematic singularity structure 
of this graph is sufficiently complex to give insight into how NDIM operates and gives some 
hope that more complicated graphs can also be evaluated. 



Our paper is organised as follows. In Sec. |27^ we first review the theoretical framework of 
one-loop integrals with Schwinger parameters and briefly explain the basic idea of integrating 
in negative dimensions. We then apply NDIM to construct template solutions for arbitrary 
one-loop integrals together with a linear system of constraints that relates the powers of the 
propagators in the loop integral to the summation variables. The system of constraints has 
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many solutions and each one must be inserted into the template solution, yielding a sum 
over fewer variables that can be identified as a generalised hypergeometric function. This 
method gives simultaneously all the solutions in all the possible kinematic regions. The 
approach is illustrated for the massive bubble integral where we show how to recover the 
known results. We discuss how the form of the solution in different kinematical regimes is 
dictated by the convergence properties of the hypergeometric functions and the structure 
of the system. In Sec. ^ we consider one-loop triangles and give the form of the template 
solution and the system of constraints with arbitrary powers of the propagators, internal 
masses, external legs off-shell and for general D. We apply this result to triangle integrals 
with three scales and give expressions valid in the various kinematic regions appropriate to 
the vertex integral. Results are given in the form of hypergeometric functions of one and two 
variables, which are defined in Appendix [A[ For specific choices of D and the propagator 
powers, these functions can be evaluated as logarithms and dilogarithms using the integral 
representations that are also provided in the appendix together with a list of hypergeometric 
identities that often simplify the results. Finally, our method is summarized in Sec. f|. 



2 Theoretical framework 



The generic n-point one-loop integral in D- dimensional Minkowski space with loop mo- 
mentum k is given by 

/„ D (Wi {Qj>, {M?» = / ^nj^r^ (2-D 

where, as indicated in Fig. [TJ, the external momenta hi are all incoming so that Yh=i &f — 
and the propagators have the form 

A 1 = 

A = Ue + kjj — M'f + iO i ^ 1, (2.2) 

Mi being the mass of the ith propagator. The external momentum scales are indicated with 
{Qi}. For standard integrals, the powers Ui to which each propagator is raised are usually 
unity. However, we wish, where possible, to leave the powers as general as possible. As 
discussed earlier, this may have some advantages in evaluating two-loop integrals where often 
one-loop integrals with arbitrary powers can be inserted into the second loop integration. 

To evaluate this integral, we introduce a Schwinger parameter Xj for each propagator 
(noting that A4 < after Wick rotation to Euclidean space) so that 




a? r(ui 



dXiXi' exp(xjAj), (2.3) 
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Figure 1: The generic one- loop graph 
and we can rewrite Eq. (|2.1| ) as 

% (W; {Qil {Ml}) =jVxJ- FJ - 2 e W (E^J > ( 2 - 4 ) 



where we have used the shorthand 



/ n 1 r°° \ 



with 



<7 = $>i. ( 2 ' 6 ) 

i=l 

The Gaussian integral over the loop momentum can be solved in a straightforward way, and 
using the Minkowski space relation 

/ ^ ex P( aA;2 ) = ^' ( 2 - 7 ) 
we have the usual Minkowski space result 

In (M; {0?}, {M?}) = J Vx^p exp(Q/V) exp(-M). (2.8) 
The quantities V and M. are given by 



V = £2*, (2.9) 
A-f = J^XiMl (2.10) 



i=i 
6 



Vi 




v 2 



Figure 2: The one- loop vacuum graph 

while Q may be simply read off from the Feynman diagram 

n-l n /j-i \ 2 9 

Q = E E ^ £aJ =E2. (2.ii) 

j=l j=i+l \k=i ) 1=1 

Each of the g terms in <2 is indicated with Qi, and is obtained by cutting the loop diagram 
into two across propagators a and b and constructing the four-momentum Qf on each side 
of the cut: Qj = x a XbQ 2 . For example, the one-loop bubble graph shown in Fig. || has two 
propagators (n = 2), so that V = x\ + x 2 and = iiMj + x 2 M|. Q is obtained by 
examining the momentum flowing across the only possible cut (q = 1): Q — Qi — XiX%Q\, 
with Qj = k\. 



2.1 The negative-dimension approach 

The crucial point in the negative- dimension approach is that the Gaussian integral ( |2.7| ) is 
an analytic function of the space-time dimension. Hence it is possible to consider D < and 
to make the definition 

^- 2 (kr = n\S n+ », (2.12) 

for positive values of n. We see that by expanding the exponential in ( |2.7|) and inserting the 
definition ( |2.12| ), after the exchange of the integration with respect to the summation 

" eMak 2 ) = £ ^ / ^- 2 (k 2 T = (2-13) 



we recover the original result, provided that D is both negative and even (so that the 
Kronecker 5 can be satisfied and the contribution with n = —D/2 selected from the sum). 
We note that, with this definition, negative- dimensional integrals can be shown to obey the 



necessary translation properties [pi] 
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2.2 The general case: different masses 



For the one-loop integrals we are interested in here, we follow the approach suggested by 
Suzuki and Schmidt [P3fl-|PE| and view Eqs. ( |2.4j ) and ( j2.8| ) as existing in negative dimensions. 
Making the same series expansion of the exponential as above, Eq. Q2.4|) becomes 



d D k " {xiAi) n * 



ni,...,n,j=0 i=l *' 



= [Vx E (-ni,---,-n»;W?},W})lI^T. ( 2M ) 

J m,...,n n =0 i=l 

where the are positive integers. The target loop integral is an infinite sum of (integrals 
over the Schwinger parameters of) loop integrals with negative powers of the propagators. 



Likewise, we expand the exponentials in Eq. 

i D n (M; {Q% Wn) = / vx E E L ^ L , (2-15) 



and introduce the integers qi, ■ ■ ■ ,q q , Pi, ■ ■ ■ ,p n an d mi, ... , m n to make multinomial expan- 
sions of Q, V and M. respectively 

oo 

Q n = E 



D 



or 

=o 9i ! 


...^V(gi + 


- + 


../..< 1 

=0 Pi ! 


rppTl 

...^r(pi + . 

pJ 


■■+Pn)! 


(- 






71=0 


mi! 


mj 


n 

Ep 

i=l 


i- -n- y 


and 



V- n ~~ = E -^...^-(pi + .-.+Pn)! (2.16) 

P!,...,p n =0 V\- Pn- 

°° C_ Tl M 2 )" 11 (-r M 2 ) mn 
{-M) m = E ( "J 7 (mx + .-. + mn)!, 

mi, 

subject to the constraints 

E^ = n ' E^ = _n_: 7 and E m i = m - ( 2 - 17 ) 

i=l i=l i=l 

Altogether, Eqs. ( |2.15| ) and ( |2.16|) give 



^ x 2^ i ii i i i \Pl ~r ■ ■ ■ ~r Pn)-, 

Pl ,.t^„=o Pi!---PJ ™i' "■ 



m„! 



, . . . ; m n — 



(2.18) 



with the constraints expressed by Eq. (|2.17|) . 
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We recall that each of the Qj is a bilinear in the Schwinger parameters, so that the target 
loop integral is now an infinite sum of powers of the scales of the process (with each of 
the Mf and the Q\ raised to a different summation variable) integrated over the Schwinger 
parameters. 

Equations ( [2.14 ) and Q2.18 ) are two different expressions for the same quantity: I® . 
However, rather than performing the integrals over the Xj's, we use the fact that the Xi s are 
independent parameters, so that the integrands themselves must be equivalent: 

oo n rii 

E 4 D (-m,...,-n n ;{g i 2 } ) {Mn)ll|7 = 

ni,...,n n =0 i=l i' 

00 Q? Q Qq x? 1 x Pn (-xiM?) mi (-x M 2 ) mn 

_ -Zl ■ ■ ■ -Z q X 1 . . .X n XjM±) ^ { X n M ) {pi+ ^ + ), (21Q) 

Pl ,.t^=o Qv---Q q ] - Pi!---Pn! rn n \ 

qi,...,q q =0 
mi ,.. .,77177,— 

Comparing term by term the left-hand-side (lhs) with the right-hand-side (rhs) allows us 
to read off the value of I®. In fact, the coefficient of the term x± Vl . . -x~ Vn in the lhs of 
Eq. (|2.19|) , where the v\ are negative integers, is given by 

^'''^(Hra^))' (2 ' 20) 

This term is equal to the coefficient of the term x~[ Vx . . . x~ Un in the rhs of Eq. ( |2.19| ). Writing 
a general expression is not possible, since the Qi are process dependent. Nevertheless, we 
can extract the momentum scale Qf from each of the Qi and find the coefficient of this term 
to be 



E (Ql) V ...(Qiy<(-M?) m \..(-M$ 



PX,...,pn— 
q\,...,q q =Q 
m-i ,...,m n =0 



1 \ (X 1 

x 



n r (1 + m ,)r(i + P,) J IS WW)) r I 1 + £N ' (2 ' 2i; 



subject to the n (process-dependent) constraints that ensure that the power of Xi on the lhs 
(— Vi) is equal to the power of x« on the rhs, which is generally a combination of integers of 
the summation. 

By adding together the first two expressions in Eq. ( [2.17| ), we obtain an additional con- 
straint, that is 

D 

Pl + ...+p n + qi + ... + q g = --. (2.22) 
Equating Eqs. ( |2.20| ) and ( |2.21| ), we obtain an expression for the loop integral with 
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negative powers of the propagators in negative dimensions 

oo 

/n D (W;{Q?},{M|}) = £ {Ql) qi ---{QT {-Mf 



mi 




Equation (|2.23|) , together with the constraints, is the main result of this paper. The loop 
integral is written directly as an infinite sum. Given that Q can be read off directly from the 
Feynman graph, so can the precise form of Eq. ( |2.23|) as well as the system of constraints. 



Of course, strictly speaking we have assumed that both Vi and D/2 are negative integers 
and we must be careful in interpreting this result in the physically interesting domain where 
the Vi and D are all positive. However, this is relatively straightforward and in the following 
sections we show how quite general results for one-loop massive bubbles and triangles can 
be obtained. 



Example: to give an explicit example of how Eq. ( [2.230 and the system of constraints 
appear, we consider the one-loop bubble with different masses. Equations Q2.14| ) 
and ( |2.18D become 

I^u^u^QlMlMl) 



Vx £ I^(- ni ,-n 2] QlMlMl)^ 



J/1 J^2 



ni,ri2=0 
oo 

Vx J2 

Pl,p 2 ,<2l,mi,ni2=0 



ni! ri2\ 

(x 1 x 2 Q 2 1 ) qi xf x p 2 2 (-xijgT (-^Mfp 
gi! pi\ p 2 \ mil m 2 ! 



(Pi +P2V; 

(2.24) 



so that, by selecting powers of x 1 Vl and x 2 1/2 , we find (see Eq. (|2.23| )) 



x 



E («?)'" {- M T (- M l) 

Pl,P2,9l,™-l,™-2=0 

r(i-u 1 )r(i-iy 2 )T(i+ Pl + P2 ) 



(2.25) 



r(i + mi)r(i + m 2 )r(i + Pl )r(i + P2 )r(i + qi ) 

together with the system of constraints 

gi+pi+mi = -v x , 

qi+P2 + m 2 = —v 2 , (2.26) 
D 

qi+Pi+P2 = - — ■ 



10 



In Sec. |2]3], we will show how this particular system can be solved to give results for 
the bubble integral in positive dimensions D, with arbitrary positive powers of the 
propagators. 



2.3 The general form of the solutions 



In general, for an n-point one-loop integral with q external momentum scales and m mass 
scales, there will be (n + q + m) summation variables and (n + 1) constraints. Altogether we 
expect (n + q+ m) \f(n+ 1) \/{q + m — 1) ! possible solutions (some of which will be eliminated 
by the specific form of the system of constraints). It is easy to see that these solutions span 
physically different kinematic regions (depending on the powers of the kinematic scales) and 
the summations will only converge in the appropriate kinematic domain. We expect that 
solutions in one kinematic region should be analytically linked to those in other domains. 

Each solution of the system of constraints, once inserted into the template of Eq. ( |2.23 ), 
has the following generic form 

TUB x SUM, (2.27) 
where we have introduced the following notation: 



SUM is the sum over the terms that contain unconstrained indices of summation. 
Instead of dealing with V functions, we have formed Pochhammer symbols, defined as 

M- 3 ^. (2-28) 

because they are the most suitable way to write generalized hypergeometric functions. 
For example, in the case where there is only one remaining summation variable n, then 
SUM takes the form 

SUM^± p4^p^L ^ (2.29) 

where x is the ratio of kinematic scales. The variables and 6j are linear in the V{ 
and D and do not depend on the summation variables. To put SUM in this form, it 
is often convenient to use the identity (see Eq. ( |2.32| )) 

fe -' !) = ( - 1) "(T^)' (2 ' 30) 

In most cases, SUM can be directly identified as a generalized hypergeometric func- 
tion, in the region of convergence of the series. In general, these hypergeometric 
functions are analytic and may be evaluated at positive values of D and z/j. 
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- The prefactor VTZS contains all the rest of the terms that are not included in SUM.. 
More precisely, it is a product of external scales raised to fixed powers, and V functions 
that do not depend on the summation variables. These may be produced either directly 
from the particular solution of the system, or in the generation of the Pochhammer 
symbols. 

In the general case of an n-point one-loop integral with q external momentum scales 
and m mass scales, inspection of Eq. (|2.23|) dictates that we produce: 

- nT functions with argument (1 — Ui), 

- (n + m + q) factorials of the summation variables in the denominator, 

- one T function in the numerator: T (1 + J2k=iPk)- 

Applying the (n + 1) constraints leaves (m + q — 1) factorials of the remaining un- 
constrained summation variables and produces an additional (n + 1) T functions in 
the denominator and one in the numerator, as Pochhammer symbols are formed using 
Eq. ( [2.281 ). Altogether there will be (n + 2) Pochhammer symbols in SUM, while VTZS 
will be a ratio with (n + 1) T functions in both numerator and denominator. In both 
SUM and VTZS, the number of functions may be reduced if there are cancellations 
between numerator and denominator. 



For physical loop integrals with positive powers of propagators, we need to evaluate VTZS at 
positive values of the 2/j and positive D. A problem is immediately obvious: the numerator 
of VTZS contains T(l — z/j), so that, for positive integer values Ui, it appears that we need to 
evaluate the T functions for negative arguments, where they are singular. However, VTZS is 
an analytic function and these singularities cancel between the numerator and denominator. 

In fact, it can be easily shown that, starting from the identity 

T(z + 1) = zT(z) } (2.31) 

we have 



riV) , T(n + l-z) 



r(z-n) v ' r(i-*) ' 

where z is a real (or complex) number, and n is a positive integer. 



(2.32) 



In the product of V functions in the numerator and denominator of the VTZS term, we can 
make an iterated use of the identity fl2.32|) , provided we treat D/2 as an integer, as we have 
already done in the multinomial expansion. We can then rewrite the T-function prefactor in 
a more amenable way by flipping all of the T functions from numerator to denominator and 
vice versa 

T7 ^ - f-i^c*-> ff r(i-fl) (2 33) 
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where the index i runs over all (n + 1) T functions in the numerator and denominator of 
VTZS. In addition, it can be shown that 

n+l Y) 

E(A = (2.34) 

which is independent of the z/j. 



2.3.1 An example: the massless bubble 

Returning to the example of the one-loop self-energy diagram introduced in Sec. |2.2| , and 
setting the masses of the internal lines to zero, Mi = M 2 = (which is equivalent to 
terminating the series in mi and m 2 at the first term), we obtain the simpler system of 
constraints (see Eq. (|2.26| ) with mi = m 2 = 0) 



qi+P2 = -»2, (2.35) 



qi + Pi + P2 



D 
~2 ' 

Since m = 0, q = 1 and n = 2, we expect that the (n + 1) = 3 constraints exactly determine 
the (m + q + n) =3 variables. In this case SUM. = 1 and the result is entirely given by the 
prefactor V7ZS. Solving this system yields 

D 

D 
~2' 
D 

P2 = Vl-~. 

Inserting these values directly into Eq. ( |2.25| ) with mi = m 2 = we find 

i?(vi,v 2 ;Qlo,Q) = Pn£ 

r (i - ui) r (i - v 2 ) y{\ + v 1 + v 2 -d) tt-vx-v. 



qi = -^-v\- "2, 



Pi = V2~- 



r (l + v x - § ) r (i + u 2 - § ) r (1 + f - v x - v 2 



(2.36) 

As expected there are (n + l) V functions in both numerator and denominator and further- 
more the arguments satisfy Eq. ( |2.34|) . We therefore apply Eq. ( |2.33|) and find 

I 2 (vi,v 2 ,Q x ^)-(-l)* T (ui)T(u 2 )T(D-u 1 -u 2 ) W ' (2 ' 37) 

where vi and v 2 are positive and which agrees with the known result straightforwardly 
obtained using Feynman parameters. 
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2.4 Massive bubble integrals 



We want now to give a detailed description of how to build the solutions starting from the 
general form ( |2.23| ) for the loop integral and from the system of constraints, and we want 
to discuss how the solutions of the system of constraints need to be combined to give a 
meaningful result. 

We will refer to a precise example to make things clearer: the bubble integral with differ- 
ent masses in the propagators of Eq. ( |2.25| ). Particular cases with u\ — v 2 — 1 are important 



in electroweak renormalization and have been known for some time (see for example [||) 
The more general cases with ^ / i/ 2 ^ 1 have been studied by Boos and Davydychev |13| 
using the Mellin-Barnes integral representations. 



As discussed at the end of Sec. |2.2| , in the case where the two masses are non-zero and 
different, there are (m + n + q) = 5 summation variables: two for the propagator masses (mi 
and m 2 ), two for the expansion of V (pi and p 2 ) and one for the external momentum scale 

(ft)- 

The (n + 1) = 3 constraints are given in Eq. (|2.26|) . There are a maximum of 5!/3!/2! = 10 
possible solutions, one for each of the ways in which we can choose three variables among the 
five, and solve the system with respect to these triplets. In this case, there is no solution if we 
try to solve the system for {pi, qi, m 2 } or {p 2 , gi, mi}, so that we have only eight solutions. 

Each of the eight solutions corresponds to different values of the integer summation 
variables and we insert each of them into the general expression for the propagator integral, 
Eq. ( |2.25| ). For example, solving for {pi,p 2 ,Qi}, yields 



D 

pi = v 2 + m 2 - —, 
D 

p 2 = fi+mi-—, (2.38) 
D 

qi = —-v 1 -v 2 -m 1 - m 2 , 



and the contribution of this solution to the integral (|2.25|) is 

j{mi,m 2 } / 2 \f-^-^-mi-ma (_ M 2\ mi (_ M A m2 T (1 — V\) V (1 - V 2 



mi ,7712=0 



X 



T(l +mi)T(l +m 2 ) 
r (1 + v\ + v 2 - D + mi + m 2 ) 



r (l + v 2 - f + m 2 ) T (l + v 1 - f + mi) T (l + f - z/! - v 2 - mi - m 2j 

(2.39) 

where we have labelled the integral with respect to the indices of summation and we have 
dropped the functional dependence of I 2 , for ease of notation. 
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As discussed in Sec. |2.2| , we now form the Pochhammer symbols, and we make use of 
the Eq. ( |2.30| ) to flip the Pochhammer symbol in the denominator with negative indices of 
summation, to obtain 



r {m 1 ,m 2 } 



r (l - t/j) r (i - u 2 ) r (1 + v x + v 2 - d) 
r (1 + p 2 - § ) r (1 + Vl - § ) r (1 + § - Vl - v 2 



™{l + Pi+P 2 -D,mi + m 2 )(p 1 + v 2 -£im 1 + m 2 ) (Mf/Ql) mi (M|/Q? 

x 2_ 



2\ m 2 



mi ,m2=0 



1 + v 2 - m 2 ) ( 1 + V\ — y, mi 



mi! 



m 2 ! 



(2.40) 



so that we can recognize the general form of Eq. (|2.27|) : the first line of the rhs is VTZS 
while the second is SUM.. By flipping the V functions in the prefactor term VTZS, using 
Eq. (^33j) , we get 



r {mi,m 2 } 



r(f--i)r(f-, 2 )r(, 1 + , 2 -f 



r^r^rp-i/!-^) 



2~ : 



x F 4 ( 1 + vi + i/ 2 - D, v x + i/ 2 - — , 1 + ui - — , 1 + v 2 - — , I , (2.41) 



D Ml Mf 
2"' Of' Of 



where we have used the definition of Appell's F4 function given in Eq. ( |A.6|) . 
In the same way, we can obtain the other seven solutions: 



7-{pi,mi} 
J 2 



(-!)*(<«) 
x F 4 ( 1 + z/i 



-M 2 ^ 2 



-i/2 



D D D 



r{p2,ni2} 



■1)# (Q\ 



— U2 



-Ml 



f r fa " f ) 



D D 

x F 4 I 1 + v 2 - —, u 2 , 1 + — - v u 1 + v 2 



D Ml Mf\ 

2"' of' ofy ' 



r{pi.P2} 



(-D#(g?)- f (-m*)^ 1 (-m|) 



# _, 2 r fa - f) r fa - f ) r (f) 



r(z/ 1 )r(z/ 2 )r(o) 



x f 4 ( 1, — , 1 + — - ^1, 1 + — - ^2, n2 ■ n2 1 



~2 



Of' Of 



r{gi,m2} 



1)T -M^ 



V\ + v 2 



D 



r(f-, 2 ) 



rwr 
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, D D D Q\ M|\ 

X Fa Z/i + Uo , Z/ 2 , ,l + l/ 2 , K , n , 

1 I 1 ; 2 2 2 Ml Ml' 



if" 1 = M)*(-^>^) ? ~"%^ 

( D D D Q\ Mf\ 

~ ( } ( 2) i> 2 )r(f) 



„ ( D D D Ml QY 



^ D D DM? Ql\ , 

x F4U2, — ,1 + i^,— ,— 7,-^ ■ 2.42 

V 2 ' 2 ' 2 ' Mi' Mi) K ' 

Solution fj pl,P2 ^ deserves a comment. Before any flipping of the T functions between 
numerator and denominator, the prefactor VTZS, has a T (1) in the numerator, due to the 
fact that the Appell's F4 function has its first argument equal to 1. Flipping this T (1) 
according to Eq. (|2.33|) generates a T (0) in the denominator, so that this solution is to be 
considered to be equal to 0. 



2.4.1 Identification of the groups of solutions using the convergence regions 



The Appell's F 4 (a, /3, 7, 7', x, y), defined in Eq. ([A.6Q , is convergent only if (see Table [l]) 

\V*\ + \y/y\<l, (2.43) 

so that we can form three different groups, according to the kinematic region of convergence 
of the series 



T?(v 1 ,v 2 ;Ql,MlM i 



I2 (fi,^5 QlMlM, 
^{u^-QlMlM, 



T {m 1 ,m 2 } 1 T {p2,m 2 } , r{pi,mi} 
1 2 ' 1 2 ' 1 2 

j{qi,rn 2 } _|_ j{pi,qi] 
j{qi,mi} _|_ j{p2,qi} 



if 
if 
if 



Ml + J Mi < JQl 



<Q\ + J Mi < J Ml 



Ml < J Ml 



(2.44) 

We note that from the convergence properties of the F4, if it was not elim inat ed b y the zero 
in the prefactor, I^ P1 ' P2 ^ would belong to the first kinematic region, J Ml + J Mi < yQf. 
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In this way, using NDIM we have simultaneously obtained all the different forms of 
the hypergeometric functions that express the integral I 2 for different kinematic regions of 
Mf, Mf and Q\. These results for I 2 agree with those obtained with the Mellin-Barnes 



method 13 



It must be noted that we can go from one kinematic region to the other, just by applying 
formula ( |A.51| ) of the analytic continuation of the F4 function. As stated at the end of the 



previous section, the appearance of T (0) in the denominator, that occurs during the process 
of analytic continuation, just kills that term. 



2.4.2 Identification of the groups of solutions using the system of constraints 

We would like to address here a different method to form the groups of solutions. This is 
based only on considerations of the system of constraints, and more precisely on the sign 
of the summed indices of the series {pi, q^, mj, without any knowledge of the region of 
convergence of the specific series. 



For this purpose, the actual value of Ui and of D in the system ( 2.26Q is irrelevant 



because it only modifies the sign of a finite number of the summation variables of the series. 
For example, solution ( |2.38| ), obtained solving the system with respect of the two indices 
{mi, m 2 }, tells us that the "bulk" sign of pi is equal to that of m 2 , because for m 2 sufficiently 
large, the contribution of v 2 — ^ is no longer important. The same thing happens for p 2 and 
qi, whose "bulk" sign is equal to that of mi and —mi — m 2 , respectively. 

For this reason, instead of considering the full inhomogeneous system, we consider the 
homogeneous one, obtained by setting ^ = and D = 0. The system ( |2.26| ) for the massive 
bubble then becomes 



qi+pi+mi = 0, 
qi+p 2 + m 2 = 0, (2.45) 
qi+Pi+P2 = 0. 

We would like to stress the fact that the last equation has always the same form, since this 
is the constraint expressed by Eq. ( [2.22|) . 



We can now build a table of signs for pi, qi and m ; . The last equation gives rise to one 
of the following cases: 
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> 


< 


Pi, P2 


Qi 


r 1 , ^Il 




P2, Ql 


Pi 


Pi 


P2, Ql 


P2 


Pi, Qi 


Qi 


Pi, Q2 



Using the last equation of (|2.45 ) to eliminate q\ from the other two equations of the 
system ( p. 45 ), we have 



mi = P2, 
m 2 = pi, 

so that we can complete the previous table in the following way: 



(2.46) 



> 


< 


pi, p 2 , mi, m 2 


Qi 


Pi, Qi, m 2 


p 2 , mi 


P2, Qi, mi 


Pi, m 2 



where the last three lines have been neglected, being equal and opposite to the first three 
ones. 

We know that the summation indices of the series must be positive integers and we can 
therefore imagine solving the system ( f2.26j ) with respect to any pair of variables that are 
simultaneously positive. The table provides us with this information and we can directly 
read from the table which integers are simultaneously positive and use them to form a group 
of solutions with similar properties by selecting all possible pairs of summation indices from 
the list of positive indices. 

Starting from the first row of the table, and considering the identities ( |2.46|) , that embody 
the fact that we cannot solve the system with respect to the pairs of indices {pi,m 2 } and 
{p 2 ,mi}, because they are linearly dependent, we can form the following subgroups: 

Pi, P2, mi, m 2 {pi,p 2 }, {Pi,mi}, {p 2 ,m 2 }, {m 1 ,m 2 }. (2.47) 

The same thing can be done with the other two rows of the table: 

Pi, q u m 2 =^ {pi,qi}, {qi,m 2 }, (2.48) 
p 2 , q u mi {P2,Qi}, {Qi,mi}. (2.49) 
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These are exactly the groups obtained by adding solutions according to their region of 
convergence (see Eq. ( |2.44j )), once we consider the fact that I^ P1 ' P2 ^ = 0. 

This method gives the correct groups only for the cases where the homogeneous system 
can be solved without any ambiguity. There are examples, and we will meet one in Sec. |3.2.2| , 
where the sign of some indices of the series are undetermined, because they have a dependence 
on other indices of the type 

Pi = p 2 +Ps, 

where P2 and p^ have opposite sign. In this case, we cannot say if pi is positive or negative. 

Up to now, we do not have a way to deal with these cases directly from the system of 
constraints, and we leave the task of further investigating this issue to future works. 



2.4.3 The limiting case: M t / 0, M 2 = 

We conclude this section, by considering some extreme cases. First we consider the limit of 
one massless propagator in the self-energy diagram. We can compute this integral in two 
different ways. 



1. We can start with the system (|2.26|) with m 2 = 0: we have four variables and three 
constraints, so that we end up with a single-index series, that turns out to be a Gauss' 
hypergeometric 2-^1 function (see Eq. (|A.1| )). 



2. We can simply take the limit for M 2 — >0 of the general expressions ( j2.41|) and ( 2.42j ). 



We can apply this limit only to the solutions that are convergent in the new kinematic 
regions, and we cannot take the limi t for solutions J^ qi ' mi ^ and j| P2 ' 91 ^ because they 
are defined only for ^Qj + Jm] < Jm$. 

The expression for F 4 (a, /3, 7, 7', x, 0) is easily obtained from its definition ( |A.6[ ) with 
the second summation series collapsing to its first term 

F A {a, (3, 7, 7', x, 0) = 2 Fx {a, (3, 7, x) . (2.50) 

Both procedures give the same result. 
If Ml < Q\ 
I?(w,Cft,M*,0) = lt l} + lt } 



D_ Ui _ U2 



r (z/x) r (v 2 ) T{D-u 1 - u 2 ) 



( D D M 2 ~ 

X 2^111 + ^1 + ^2-^,^1 + ^2-^,1 + ^1-^,-^2" 
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+ {~ l r (Qi) [-Ml 



D D Mf 

2' 1+ 2~ _z/1 ' gf 



x 2^(1/2,1 + 1/2 + , (2.51) 



if Q? < M 2 

if (z^q^m^o) = 7; 



{51 } 
2 



(-l)f f-M^ ^ r (^ + ^-T) r (T-^ 



/ D D Q\ x 

x ^[i^ + i*--,!*,-,^]. (2.52) 



2.4.4 The limiting case: — > 



Similarly, we can take the limit of the general self-energy diagram where the external mo- 
mentum scale vanishes. Once again, we can either return to the system ( |2.26| ) with one fewer 
variable (qi = 0) or we just take the Q\ — »0 limit of the general result ( 2.44|) in the appro- 



priate kinematic regions: yQj + y Ml < y 1 Mf or yQf + y Mf < yMf. Both procedures 
yield the same result: 

If M l > M 2 

I 2 D (v 1 ,v 2 ;0,M* 1 ,Ml)=lt 2} + lt } 



1)f (_ M 2)f — r(, 1 + , 2 -f)r(f-,2 



x 2 F 1 [u 2 ,u 1 + u 2 --,l + u 2 --, J ^ 



D 



£-„r(i»,-§) / D D Mf 



+ (-l)t ( -^)"-(-^)T-'_y_ii lFl ^- |1 + -_^^J | (2 . 5 3) 
with the result for M 2 > Mi obtained by the exchanges Mi <-> M 2 and v\ <-» v 2 . 

Provided that we do not violate the validity of the kinematic regions, we can take the 
subsequent limits of the energy scales. For example, we can safely take the M 2 — > limit for 
the solution where M 2 < M\. In this case, only the first term survives and we obtain the 
familiar result 



„ d r U + u 2 - § ) r (f - u 2 ) 

I 2 > 1; v 2 - 0, Ml 0) = (-I)t (-Ml) 2 -V 2 J \ 2 i. (2.54) 

r K) r (yJ 
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2.5 The special case: all masses equal 



For the special case where each propagator has the same mass, Eq. ( [2 -ICQ becomes 

n 

M = Y, XiM 2 = VM 2 , (2.55) 



i=i 



and we have an important simplification. As before, we expand the exponentials in Eq. ( f2~8 
and make multinomial expansions of Q and V 



n=0m=0 nlm[ 



00 Of Q q i r? 1 r Pn (-M 2 ) m 

vx E , ? V , V (pi + --- + Pr. 



P1 ,."=o Qr-'-Qq'- PV----Pn<- mi 

m = 



(2.56) 



subject to the constraints 



q n D 

2_. Qi = n an d ^jPi = ~ n V m. (2-57) 

i=l i=l ^ 

Equating Eqs. ( |2.14| ) and ( |2.56| ) and, once again, identifying powers of x~ Ui , we obtain 
an expression for the loop integral with negative powers of the propagators in negative 
dimensions for all masses equal 



oo 

4 D (M;{^ 2 },M 2 ) = £ (Qi) qi ...(Ql) qq (-M'< 



Plf-,Pn=0 
m=0 



X 



l[r(i+ Pi )){l\r(i + qi )) r(i + m) ' ^ 



subject to n constraints that each of the powers of Xi match up correctly. However, the 
constraint that matches up the powers of Q and V, obtained by summing the two expressions 
in Eq. (|2.57|) , is now 

D 

pi + . . . + p n + q x + . . . + q q = -— + m, (2.59) 

rather than Eq. ( ggg ). 

We see that there are (n + q + 1) summation variables and (n + 1) constraints, leaving 
q remaining summations. We note that the structure of the solution is precisely as for the 
unequal-mass case and is treated in the same way by constructing the sum over Pochhammer 
symbols SUM. and the Y function pref actor V7ZS. 
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Example: to give an explicit example, we consider the self-energy correction to the 
propagator integral with equal masses. There are (n + q + 1) = 4 summation variables 
with (n + 1) = 3 constraints. In this case, the template solution Eq. ( j2.5S| ) is given by 



pi,P2,qi,m=0 

x r(i-^)r(i-^)r(i+ Pl +p 2 ) 
r(i+ ft )r(i+ ft )r(i + ft )r(i + m)' 1 ' ] 

while, matching the powers of Xi gives the system of constraints 

qi+Pi = -vi, 

qi+P2 = -u 2 , (2.61) 
D 

qi+Pi+P2 = -—+m. 

There are four summation variables (pi, p 2 , q\ and m) and three constraints, and we 
obtain four series solutions, with only one index of summation. 

Defining a = V\ + v 2 (see Eq. (|2.6|)), we have: 
If Qj > AM 2 

I? (u u v % Ql M 2 , M 2 ) = lt } + 4 P1} + 4 P2} 



■D* (Ql 



(7 



D 1 a D D D D AM r 



X 3^2 1 + - - — , - + - - — , O ~ — , 1 + V X - — , 1 + V 2 



+ (-i)M^ (-M 



2 2 ' 2 2 2' 2' 2' 2 ' Q? 

f r (u 2 - f ) 



„ ( v\ v 2 1 v x v 2 D 4M 2 

x sF 2 U,l + - --,- + - --,1+^-^,1 + --^,-^ 



4 ,• 1> 



v 2 v\ \ v 2 V\ D 4M 



2' 



x 3 F 2 \u 2 , 1 + — - — , - + — - — , 1 + v 2 - u u 1 + — - Vx 



2 2 '2 ' 2 2' ' A 17 ' 2 x ' Qf 



(2.62) 



if Q 2 < 4M 2 

If U,v 2 ;Ql,M\M*)=lt } 
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-M' 



2\ 2 



U1 - U2 



r (--f) 



(2.63) 



In forming the Pochhammer symbols we have made use of the following duplication 
formula 



z,2n)=4 n (-,n) - + -,n . 



(2.64) 



The procedure described in Sec. |2.4.2| on how the solutions group together to give the 
correct answer in a particular kinematic region, is straightforward. In fact, the solution 
of the homogeneous counterpart of the system of constraints ( |2.61[ 



is 



p 1 =p 2 = m = -qi. 



(2.65) 



We can then form groups for the indices of summation of the series out of combinations 
of indices with the same sign, that is: 



{Pi}, {P2}, {m} or {<?!}. 



;2.66) 



Not surprisingly these were the groups of solutions formed by considering the conver- 
gence properties in Eqs. ( [2.62j ) and (|2.63|) and reproduce the known result |fL3|j .. 



3 Massive Vertex integrals 

We now turn to massive-triangle integrals where each propagator can have a different 
mass and each external leg can be off-shell. The propagators and momenta are labelled as 
in Fig. |3|. Throughout this section, the number of propagators n is equal to three and in the 
most general case, we have 



V = xi + x 2 + x 3 
Q = x 2 x 3 Q\ + x?,xi Ql + xix 2 Ql 
M = xiMl + x 2 Ml + xzMl, 



(3.1) 



where Q\ = kf. 



Based on the discussion of Sec. |2]2|, the NDIM method provides a generic solution with 
(n + q + m) =9 summation variables and (n + 1) =4 constraints, with the template solution 
given by 



rD 



(u u u 2 ,u 3] QlQlQl,MlMlM. 
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Figure 3: The one- loop vertex diagram 



E «' («9 (<s) (- m (- m >) (-< 

r(i - f ( ) 



P1,--,P3=° 
91, ■■•,93=0 
, .. .,7713— 



1 1 r(i + m i )r(i + ft )r(i + g i ) / 

while the system of four constraints is 



r(i+ P i+p 2 + p 3 ), (3.2) 



q2 + qz + pi + mi = -vu 
qi + q3 + P2 + m 2 = -u 2 , 



qi + q2+ps + m 3 = -z/ 3 



D , , 

Pi +P2 +P3 + 9i + ?2 + <?3 = -y- (3.3) 

The number of possible solutions satisfying this system is 9!/4!/5! = 126 of which 45 are 
eliminated by the particular nature of the system leaving 81. As usual, insertion of these 



solutions into the template yields contributions of the general form ([2.271) , where SUM. is a 
product of Pochhammer symbols and ratios of energy scales summed over the five remaining 
variables. The prefactor VTZS vanishes in a further 12 instances, leaving 69 solutions which 
are distributed among the various kinematic regions. 

At present, the technology for dealing with five-fold sums (and their integral representa- 
tions) is not sufficiently developed to handle the completely general case. For the remainder 
of this section, we therefore concentrate on particular cases of the vertex integral where 
some of the energy scales vanish, leading to either single or double sums which have been 
well studied. 
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3.1 Massless propagators: M 1 = M 2 = M s = 



We first consider the special case where all of the internal lines are massless, so mi = m 2 = 
m 3 = in Eq. ( |3.3| ), leaving six summation variables. Of the 6!/4!/2! = 15 possible solutions 
of this system, three are eliminated by the system, leaving twelve. The three-mass triangle 
is an extremely symmetric system and there are three allowed phase space regions: 



region I : 
region II : 
region III : 



'Ql> 



'Ql + \M 



> 



> 



'Qi + 



!2i 



(3.4) 



with each region bounded by the external phase-space constraints 

A 3 (Ql,Ql,Ql) >0, (3.5) 

where 

A 3 (x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. (3.6) 

The twelve solutions populate the three kinematic regions equally, four in each, and are easily 
identified as belonging to a particular region by studying either the convergence properties of 
the double sum or by considering the system, as in Sec. [2.4.2 . For example, the four solutions 



belonging to region I (qi negative, q 2 and q 3 positive) are those where the summation variables 
include pairs in the set {p 2 ,p 3 , q 2 , q 3 }, that is {p 2 ,p 3 }, fe,^}, {^2,^3} and {q 2 ,q 3 }, where 
{p 2 , q 2 } and {p 3 , q 3 } have been eliminated by the system. 



As usual, each solution is inserted into Eq. ( ^.2[ ) and treated according to the procedure 
described in Sec. |2.2j : the summation variables are converted into Pochhammer symbols; the 
T-function prefactor is flipped using Eq. ( [2.33|) and the remaining summations converted into 



generalised hypergeometric functions. In each case, we identify Appell's F 4 (a, (3, 7, 7', x, y) 
function (see Eq. |A.6j ), which, according to the convergence criteria of Table [l|, is well defined 
when yfx + ^fy < 1, precisely matching on to the physically allowed phase space. 

Summing the four solutions we find, in the region \f~Ql > \J~Q% + \fQl, 

I 3 D (u,, u 2 , v 3 - Ql Ql Ql 0, 0, 0) = lt q3} + 4 P2 ' 93} + 4 P3 ' 92} + 4 P2 ' P3} 

n 



-1)t iQ 2 ^—^ 



r(f-, 1 -, 2 )r(f-, 1 -, 3 )r(a-f) 



T(u 2 )T(u 3 )T(D-a) 

_ / D D D Q 2 Q 2 \ 

x F,[v 1 ,a--A + u 1+ v 3 - - l + ^ 1 + ^ 2 -- _ _j 

V ' v ' r \v\) r (f 3 ) V [D — a) 
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D D D Ql Q\ 

!/ S ,H U-l —I/3.I + I/1 + Uo , "TTTT , 



D 



2\ 2 



^3 



r(f-, 1 -, 3 )r(, 1 + , 2 -f)r(f-, 2 



r(z/ 1 )r(z/ 2 )r( J D-a) 

/ D D D Ql Ql 

X F 4 I 1/3, — - u 2 , 1 + V! + 1/3 - y, 1 + y ~ "1 ~ "2, g| > q2 



r (^1 + ^ 2 - f ) r + ^ - f ) r (f -v x ) 



r(z/ 1 )r( I / 2 )r(z/ 3 ) 

n ( ^ D D D Ql Q\\ , „ N 

x F 4 L> - a, - - v x , 1 + - - v x - 1/3, 1 + - - v x - v 2 , ^| , (3.7) 



which agrees with that obtained by Boos and Davydychev |L3[ using the Mellin-Barnes 
integral representation. Similar results are obtained for the other two kinematic regions, 
either by directly summing the solutions valid in that region (pairs from {p x ,p 3 , q x , 93} or 
{pi,P2, Qi, 92}, with {pi,qi}, {^2,92} and {]?3, q 3 } excluded by the system) or by analytic 
continuation of the Appell's F4 function using formula ( |A.51|) . 

Note that if one of the v t vanishes (equivalent to propagator i shrinking to a point), only 
a single term remains. For example, if v\ = 0, only the first term of Eq. ( |3.7|) survives (the 
others being killed by 1/T (0)), and the Appell's function collapses to 

F 4 (0,/3,7, 7 / ,x,i/) = l, 

as can be seen from the definition ( |A.6|) , yielding 

J 3 D (0, v 2 , u 3 ; Ql, Ql, Ql 0, 0, 0) = J 2 D (v 2 , u 3 ; Ql 0, 0) , (3.8) 

as it should. 

We can obtain some other interesting limits if we set to zero one or two external invariants. 



1. One light-like external momentum: if the zth external leg is light-like (Qf = 0), 

we can return to the general case and solve the system (|3.3| ) with = 0, or we can take 
the appropriate limit of the general solution. These limits can be safely made provided 
that we start from a valid kinematic region. In the region of validity of Eq. ( |3.7|) , that 
is \fQl > \J~Ql + \fQl, we can surely take the limits for Ql — > or Q\ — > 0. In this 
last case, for example, the last two terms in Eq. (|3.7|) vanish, while the first two terms 
collapse to Gaussian hypergeometric functions, according to Eq. ( |A.20| ), yielding, in 
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the region Q\> Q% 

iffa, u 2 , v 3 - Ql Ql 0, 0, 0, 0) = J 3 {<?2} + I 3 {p2} 



l 3 

r (f - Vx - V 2 ) T(f -V X - V 3 )T(v X + 1 y 2 + U 3 -^) 



+ (-l)T (Qi) (Q 



T(u 2 )T(u 3 )T(D -vi-u 2 - z/ 3 ) 

( D D Ql\ 

x^^a- -1 + ^ + ^3-- _j 

r (f - m - z/ 2 ) r (f - z/ 3 ) r + ^ - f 



t(i/ 1 )t(u 3 )t(d - ^ - z/ 2 - v 3 ) 



D D Q 



x 2 F X \u 2 , - - z/ 3 , 1 + — - u x - p 3 , g|J . (3.9) 

Analogous results valid in the region Q\ > Q\ can be obtained either by starting 
from the expression for I® in region II , that is \J~Q% > \J~Q\ + \J~Ql, or via analytic 
continuation of Eq. (|3.9|) , according to Eq. ( |A.49| ). 

2. Two light-like external momenta: in a similar way, we can obtain the result for 
two light-like external momenta, Q 3 = Q\ = for example, by simultaneously taking 
both Ql and Q\— >0 in Eq. (|3.7|). Only the first term in Eq. Q3.7| ) survives, and the 



Appell's function collapses to 



F 4 («,/?, 7 , 7 ',0,0) = 1, (3.10) 



D 
2\ 2 ' 



yielding 

if (y^v^Ql, 0,0,0,0,0) = (-l)S(Q? 

x r (f-^-^) r (f-^-3)r(^-f) 
r (z/ 2 ) r (v 3 ) r p - a) ' 

which again agrees with the known result straightforwardly obtained using Feynman 
parameters. Alternatively, we could have returned to the general system ( |3.3| ), where, 
with only the Q\ scale (g 2 = qz = m i = m 2 = ITI3 = 0), we would have had (n+q+m) = 
4 summation variables and (n + 1) =4 constraints. 



3.2 Two massive propagators and one off-shell leg 

We now turn to triangle integrals with two internal mass scales and one external scale. These 
have (n + q + m) = 6 summation variables and (n + 1) =4 constraints and are therefore 
described by double sums. 
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3.2.1 Mx = 0, Ql = Ql = 



In this case, the system of constraints is obtained by setting m\ = q 2 = q 3 = in Eq. (|3T3|). 
The first constraint is simply p\ = —v\ and there are only 8 solutions for the system. As 
usual, each solution is inserted into Eq. (|3.2|) and treated accordingly to the procedure of 



Sec 2.2. 



The solutions can be grouped either by studying the physical thresholds of the integral (or 



the convergence properties of the series) or by considering the system, as in Sec. [2.4.2 
The threshold for the production of two massive propagators on-shell, \JQ\ 



M 2 + M 3 , 

becomes evident when we inspect the hypergeometric functions: four solutions are convergent 
'Ml + jQj and 



above threshold, 



Ml 



M| , while the other solutions equally populate the regions 



'Mf> 



Mf > 



'Ql 



Consideration of the system reveals that the first group of solutions are pairs from the set 
{m 2 ,m 3 ,p 2 ,p 3 } and the other groups formed are from the sets {qi,P3,m 2 } and {qi,p 2 ,m 3 } 
respectively. The apparent overlap between the groups, solutions formed from the pairs 
{m 2 ,p 3 } and {m 3 ,p 2 } are excluded by the system. In each case, we identify Appell's F 4 
function, whose convergence properties match onto the anticipated regions. 

We find: 



If > y/Mf + yjM$ 
I°{p u v 2 ,p 3 ,Ql 0,0,0, Ml M { 



r {m 2 ,m 3 } 



{m 2 ,p 2 } | r{m3,P3} 
t J- 3 



+ 1 



{P2,P3> 



:-i)- (Qi 1 



D 



- T -v 1 -V2-v i 



r (y x + v 2 + v 3 - § ) T (f - ^ - v 2 ) Y(^-v x -v 3 
r (v 2 ) r (u 3 ) Y (D — Vi — v 2 — v 3 ) 

D D Ml Mf 




+ (-i)^ [Qi 



-M. 



D 



,r u+v 3 -% r £ -v x -v 2 



r> 3 )r(f -v 2 



D 



D 



m? Mr 



+ 



x F 4 \l + v 2 - —, v 2 , 1 + v x + v 2 - —, 1 + - - v x - v 3 , 

^§-»x-» r U + v 2 - f ) r (f - ut - v 3 , 



(Qiy Ki (-Ml) 



r> 2 )r(f -v 3 




+ (-1)' (Qi 
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X 



X 



r(z/ 1 )r(z/ 2 )r(z/ 3 ) 



D D 
F 4 I 1 - vt, — - v x ,± + — 



D 
~2 



M? .Ut 



(3.12) 



if v/mi > y/CR + y/M2 

tffa, V2, i*| Ql 0, 0, 0, Ml M 2 ) = 7 3 { ™ 2 ' 9l} + 4^ 



= (-l)-(-M 3 2 ) 



f-^-^-^ r (u x + u 2 + v 3 - § ) r (f - z/i - z/ 2 



r> 3 )r( 



-1^ 



-Mi 



2 N 2—^-1/2 



(-M< 



-^3 



r(, 1 + , 2 -f)r(f-, 1 



rh)r f 



~ ( D D D Ml Q 2 \ 



(3.13) 



while the result for \j Ml > y M| + yQ? is obtained by the exchanges M 2 
in Eq. ( |£T5D . 

We can check that these expressions are valid in certain limits. 



M 3 , ^ 2 z/ 3 



Checks 



The i/i -> limit: 1° (0, i/ 2 , v 3 \ Q 2 , 0, 0, 0, Mf , 

Pinching out the first propagator, the first three terms in (|3.12|) and both terms in 
( |3.13| ) survive, yielding the general bubble integral of Eq. ( |2.44| ), in the respective 
kinematic regions, 



/ 3 D (o,z/ 2 ,z/ 3 ;Q 2 ,0,0,0,M 2 ,M 



tD 



(^ 2 ,z/ 3 ;Q 2 ,M 2 2 ,M 



(3.14) 



The i/ 2 -> limit: If 0, v 3 ; Q\, 0, 0, 0, M 2 2 , 
This limit should produce a self-energy integral with no external momentum and a 



single internal mass M 3 . This is indeed the case: only the second term in ( |3.12j ) and 
the first term of ( |3.13|) survive, each yielding the same result of Eq. ( |2.54|) . 



The M 2 — ► limit: 1° (i>-i, u 2 , is 3 ; Q\, 0, 0, 0, 0, M. 

Here the real production threshold occurs at Q\ = Mf and in this limit, Eq. ( |3.12| ) 
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provides the Q\ > Mf result, while Eq. ( |3.13| ) gives the expression for Mf > Q\. 
The Appell functions again collapse to form Gaussian hypergeometric functions (see 
Eq. (pOTD ), and we find: 



If JQI > J Mi 



I?(vi,v 2 , v 3 ;Q'i,0,0,0,0,M-£ 



/j m 3} + j{ps} 



(-i)-(Q?) 

x 2 Fi { 1 + a 



-^■-v\-V2-v-i 



To 



z^i - ^2 



v\ - v?, 



D 
~2 



r (z/ 2 ) r (z/ 3 ) r (d - a) 

D M, 2 \ 



A CT - — > 1 + v \ + ^3 - TT, tJ" 



x 2 Fi ^2, 1 



r(^ + ^-f)r(f 



^1 - V 2 



D 



D 
~2 



r (u 3 ) r (f - i*) 

Mf 



^1 - ^3, "Tvf 



(3.15) 



if y/M$ > VQl 

if fa, u 2 ,u 3 ; Ql 0,0, 0,0, Mi 



r{9l} 



(-l) f (-M 3 2 ) 



-~-V 1 ~V2-Vz 



r(.-f)r(f-, 1 -, 2 ) 
r fa) r (f ) 



2 Fi u 2 , a 



D D Ql 

2 ' 2 ' Ml 



(3.16) 



This latter result agrees with that obtained by taking the limit Q\ — > 0, Q\ — > in the 
general result given by Boos and Davydychev |13[ for a triangle loop integral with a 
single massive propagator. 



3.2.2 M 3 = 0, Q2 = Q2 = o 

We now consider the triangle graph where M 3 = and Q\ = Q\ = 0. Although this graph 
is not usually present in Standard Model processes, the analysis of this graph turns out to 
be rather more subtle than the preceding triangle integrals and we will therefore describe it 
in more detail. Inspection of the singularities present in the loop integral via the Landau 
equations reveals that threshold singularities occur at Mf = Q\ + Mf. We expect that this 
equality will provide the appropriate boundaries of regions of convergence when considering 
the convergence properties of the generalised hypergeometric functions. Furthermore, since 
the convergence properties of these functions only depend on the absolute value of ratios of 
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2 

M?/Q? 



3 



4 



Figure 4: The kinematic regions for the one-loop triangle with Q\ = Q\ = Mf = 0. The 
solid line shows the threshold in the Landau surface at Mf = Q\ + Mf, together with the 
reflections Mf+Mf = Q\ and Mf = Mf+Qj. The reflections are relevant for the convergence 
properties of the hypergeometric functions which only involve the absolute values of ratios 
of the scales. The dashed lines show the boundaries Mf = Mf and Mf = Q\. 

scales, we expect that the reflections, Mf + Mf = Q\ and Mf = Mf + Q\, will also form 
boundaries in the large Q\ and Mf regions, respectively. We also expect that, in certain 
limits, the solutions match onto the kinematic regions relevant for simpler integrals. For 
example, in the limit Mi — > 0, the discussion of the previous section informs us that the 
solutions divide according to whether or not Q\ > Mf. Similarly as Q\ — * 0, there should be 
a threshold at Mi = M 2 . 

In anticipation, we therefore divide the kinematic regions up as follows: 



region I : 



Mf > Q\ + Mf, 



region 11(a) : 
region 11(b) : 
region 111(a) 
region 111(b) 



Ql > Mf + Mf and Mf > Mf, 
Q\ > Mf + Mf and Mf > Mf, 
Mf >Q\ + Mf and Q\ > Mf, 
Mf >Ql + Mf and Mf > Q\, 



(3.17) 



as shown in Fig. [|. 
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With this set of scales, the system is given by Eq. ( |3.3| ) with m 3 = q 2 = q 3 = 0. As 
usual, we construct the solutions by solving the system, inserting the solutions into Eq. ( |3.2| ) 
and following the procedure outlined in Sec |2.2j . Labelling each solution by the summation 
variables and using the definitions of the hypergeometric functions of Sec. |A.1|, we find: 

« - r (z/x + v 2 + v 3 - f ) r (f - v x - z/ 3 ) 



/ D D D M\ QV 



r fa) r (f 



^ ( D D D M? QV 



H mi ,m 3 } .^f / o2 \f-l-^-3 r + Z/ 2 + Z/ 3 - f ) r (f - ^ - Z/ 2 ) r (f - //! - Z/ 3 

3 1 j ^ rHr^rp-.!-^-^) 



, D D DM} Mf\ 

x Si ( 1 + a - D, a - —, v u 1 + v x + u 2 - — , 1 + V\ + u 3 - — , — gj, ) , 



r{m 2 ,p 2 ] _ , ^D( n2 \~^( 



J™, = (-l)f Q} ~[-Mi 



r 


+ ^ - f ) 


r (z/ 3 - z/ 2 ) r (f - us) 


r (z/0 r (z/ 3 ) r (f - z/ 2 ) 





n / D D -D M 2 M|\ 

x 6i I 1 + z/ 2 - — , u 2 , — - u 3 , 1 + v 2 - u 3 , 1 + — - v x - u 3 , --qTi -Q2 J ' 



r (z/O r (z/ 2 ) r (f - z/ 3 ) 



, D D D Ml Ml 

x S 2 [u 2 - u 3 , u x + z/ 2 - — , 1 + u 3 - — , z/ 3 , 1 + z/ 2 - — , j^, j . 



it' ps} = (-1)* (o?n \-Miy\-Mi) 



i , r o\ _, 'i / , r .)\T-"2 r {v 2 2 



„ / D D D Ml Ml\ 

x F 3 \v x ,l + v 3 - - -u 3 ,u 3 ,l + _j, 

- - r (z/! + ^ 2 - f ) r (f - i/x - u 3 



r (z/ 2 ) r (f - z/ 3 ) 



( D D D Ml Mf\ 

X if 2 ^Z/ 1 + Z/ 2 --,Z/ 1 ,l + Z/3--,Z/3,l + Z/ 1 + Z/3- - ^, , 
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r{P2,P3> 



;-i)S(q?P(-jk?) 



#-.1-^3 



-M 2 2 



r (i/j + z/ 3 



r (1^2 - ^3) 



TT ( D D ■ D M l M l 

X # 2 [ V 2 - I/3, — - V 3 , 1 + 1/3 - — , U 3 , 1 + — - U X - V3, J^, - 



Ql 



2 / ' 



r{m 2 ,pi} 



(-!)*(«)' 

X 5 2 Z/ 2 + Z/ 3 



(-M 2 )~ 



1) 



r (iaj) r (z/ 3 ) v (d — u 2 — u 3 ) 



D D D Ml Q\\ 

-l + u 2 + Us-D,u 1 ,--Us,l + u 2 --, ww ), 



r{m 2 ,qi} 



. D 

•IT 



-My 



2 \ -j" -^1-^2-^3 



r ^1 + z/ 2 + 



^3 



D 



r 



Z7 2 - Z/3 



r (u,) r 



D D D Qj Ml 



x Si I i/ 2 , v x + u 2 + us ~ y , ^3, 1 + ^2 + ^3 - 2 , 2 , 



r{Pi-9i} 



x H 2 \u 2 + u 3 



M^ Ml)' 

r(, 2 + , 3 -f)r(f-, 3 ) 



r> 2 )r 



~2 



Mf 



(3.18) 



t2 Mi y 

We now need to study the convergence properties of these solutions. For example, by in- 
specting Table [l] in Appendix |A.1| , we see that the function S\{. . . ,x, y) is convergent when 
\ x \ + I V I < !• This implies that solution /,j mi ' m2 } j s convergent when 



-Ml 



Ql 



+ 



< 1, 



(3.19) 



or, in other words, 

M 2 + Ml < Q\, (3.20) 

independently of whether Mi is larger than M 2 or not. This series therefore converges in 
both regions 11(a) and 11(b). 
On the other hand, J^ mi ' P3 ^ converges when 



Ml 



1 1 3 

2 



M| 



+ 



Ql 



M| 



> 1 



and 



Ml 



M| 



< 1 



and 



Ml 



Ql 



< 1, 



or, alternatively, 



Ml + Ml < Ql and M 2 2 > M 2 . 



(3.21) 
(3.22) 



which corresponds to region 11(b) only. 



Applying the convergence criteria to each of the eleven solutions, we find that they are 
distributed as follows: 



in region I 
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I°(v u v 2 ,v 3 ;Ql,0,0,MlM 2 2 ,0) = jj"™> + J***>, (3 23) 

in region 11(a) 

I^(iyi,u 2 ,u 3 ;Ql,0,0,MlMlO) = 1^'™^ + 4 m2,m} + lt 2 ' m} + (3-24) 
in region 11(b) 

0,0,M 2 ,M 2 2 ,0) = J 3 ^^} + 4- 2 .M + / W,P3} + / {« >P 3} ) (3 .25) 

in region 111(a) 
in region 111(b) 

/3 D (^i^2,^Q 2 ,0,0,M 2 ,M 2 ,0) = 4-^ + 4^}. (3.27) 

We see that in region 11(a), two of the solutions ^j^ m2 ' P2 ^ an d 4 m2 ' P3 ^ contain dangerous 
T functions when v 2 = v 3 . These divergences usually indicate the region of a logarithmic 
analytic continuation and can be regulated by letting v 2 = z/ 3 + 5, canceling the divergence, 
and then setting 5— > 0. Similarly, the two divergent contributions in region 11(b) ^4 m2,P2 ^ 

and 4 P2 ' P ^) also cancel in this limit. 

We have performed several checks of the correctness of this assignment into groups. 



Checks 



Analytic continuation 

Applying the analytic continuation formulae given in Appendix [A], we can see that 
the solutions are connected to each other. For example, applying Eq. ( |A54j ) to the 
F 2 functions in region I produces the H 2 and Si functions of region 11(b). Similarly, 
Eqs. ( |A.56| ) and ( |A.53| ) transform the Si and H 2 solutions of region 111(b) into the F 2 



functions of region I. 

The i/i = v 2 = u 3 = 1 limit: If (1, 1, 1; Q\, 0, 0, M 2 , M 2 2 , 0) 

All the groups give the correct answer when all the propagators are set equal to one. 

As an example, we consider region 11(b), so that we can explicitly show the cancellation 
of the 5 poles. We fix vi = v 3 — 1, v 2 — 1 + 5 and D = 4 — 2e. For these choices 
of the parameters, the hypergeometric functions simplify using the identities given in 
Sec. |A.3.2| , and we find 



/ mim9l iV e r(i-e) 2 / Q 2 \ / M? \ 



T {m 2>P2 } N e T(l-e) T{l-S) {Qlt / 2 2 M 2\- e ~ d ( o 2q] 
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T {mi,p 3 } 
J 3 



r{P2,P3} 



Mf 



( Ql ' 

e 2 Ml -Ml \-Ml 
N e (-1)- 5 I 



F 2 

Ql 



3 e 5 (M| — Mi f \-Mf 

where we have defined 



1,1,6,6 + 1,1 



lF X 1,6,1 



M 2 



M 2 - Mf ' Q 2 
M 2 \ 



N f = T(l 



-Ql 



-1-e 



(3.30) 
(3.31) 

(3.32) 



Using Eq. (|A.5Q ), we can rewrite !■ 



{P2,P3} 



r{P2,P3> 



Ne 
c2 



Ql 

-Ml 



Ql 



in the following way 

i Qj + Ml-Ml 



Mn 



Ml 



Ml 



Mn 



N e r(l + e + 5)(-l)- d r(l-5) (Qj 



Ql + Ml 



Ml 



-e-6 



(3.33) 



(5 + e) r(l + e) 5 (-M 2 ) e 

where we have safely put 5 = in the first line, since this is a finite quantity in 5. 

We see that the poles in 5 clearly cancel between /j m2,P2 ^ an d I^ P2 ' P3 ^ (see Eqs. ( |3.29|) 
and (|3.33|) ), leaving a finite remainder that is straightforwardly obtained by Taylor 
expansion about 5 = 0. 



Up to now, we have not required e to be small and expressions (|3.28|) - (|3.33|) are valid 
in arbitrary dimension D. 

If we make the usual expansion for e — > we recover the result 

'Ql + Ml 



^(l,l,l;Qf,0,0,M 1 2 ,M 2 ,0)=iV e 



Mf 



Mf 
Ml 



u, i •-■ : \ Li 2 ( i 1 + o (e) , 

(3.34) 

where we have used a series expansion for the integral representation of the functions 
2-Fi and F 2 , given in Eqs. ( |A.26| ) and ( |A.29| ), and where Ql~^Ql + i0, to recover 
the correct prescription in the Feynman integrals. We describe the details of the e 
expansion in Appendix [A.2.1| . Expression ( |3.34j ) is finite in e, as it should be, having 
no soft or collinear singularities, despite the fact that the individual contributions 
contain poles in e. 



The i/i -> limit: If (0, i/ 2 , f 3 ; Ql, 0, 0, M\, Mf , 0J 

If we set v\ — >0, we produce a one-mass (M 2 ) bubble integral with external scale Q\ 
and internal propagators raised to the powers v 2 and u 3 . 

In the different regions ( |3.17| ) and for the different groups of Eqs. ( |3.23| )-( |3.27| ), we 
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have 



I : Ml > Ql and M 2 > M 2 

11(a) : Ql > Ml > Ml 

11(b) : Q\ > Ml > Ml 

III (a) : Ml >Ql> Ml 

111(b) : Ml > Ml > Ql 



I D 



r{mi,qi} 



z/ 1= 



i/i=0 



i/i=0 



vi=0 



i/i=0 



i/i=0 



r{mi,ra 2 } 



r{mi,rri2} 



r{Pl.P3> 



1/1=0 



1/1=0 



i/i=0 



+ I 



{mi,p3> 



1/1=0 



r {m 2 ,pi} 



1/1=0 



+ /■ 



{P1.P3> 



1/1=0 



r{pi,9l} 



1/1=0 



where we have used the shorthand notation 

J 3 D = /f (0, z/ 2 , z/ 3 ; Ql o, 0, M 2 , M|, o' 

V\ =U V > 

and where the missing terms have been killed by the T (0) in the denominator. 



(3.35) 
(3.36) 



It is straightforward to evaluate the different solutions when v\ = 0. In fact, taking 
j{mi,qi} ag exam p] e; we can use reduction formula 



F 2 (a, 0, /3', 7, 7, x, y) = 2 Fi (a, /3', 7', y) 



to recover 



/ 2 D (^,^;Qi,M 2 2 ,o) 



(3.37) 



(3.38) 



in region I, that is Eq. ( p.52| ). The same thing happens to the solution in region 111(b). 

The other part of the bubble integral valid when Ql > Ml is produced by the other 
solutions in the region 11(a), 11(b) and 111(a), and agrees with Eq. ( |2.51|) . 



The v 3 -> limit: 1° (u u u 2 , 0; Ql, 0, 0, M\, M|, 0J 

Likewise, we can set z/ 3 — > producing a two-mass bubble (Mi and M 2 ) with external 
scale Ql = 0, for which the result is given in Eq. fl2.53|) . We could repeat the reasoning 
made for the previous case, and build a table of surviving solutions, analogous to ( 3.35Q . 

For example, the two terms in Eq. ( |3.23|) collapse to form the correct Gauss' hy- 
pergeometric functions when M| > M 2 and M| > Q\. Similarly, the result when 
Ql > Ml > Ml is produced by the third and fourth term of Eq. ( |3.25| ), for region 
11(b). 



- The Mx -> limit: I" \u x , u 2 , u 3 ; Q{, 0, 0, 0, M 2 2 , 

Here we expect to reproduce the result for the triangle integral given in Eqs. ( |3.15|) 
and (|3.16|) , with the exchanges M 3 <-> M 2 and z/ 3 <-> z/ 2 . Clearly in regions 11(a), 111(a) 
and IH(b), it is inappropriate to take this limit, since M 2 > Mf . In fact, if we just go 
ahead and apply the limit blindly to the solutions for regions 111(a) and 111(b), we just 
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obtain zero. 

On the other hand, in regions I and 11(b) it does make sense to send Mi — > since Mi 
is allowed to be the smallest scale present: 



I : Ml > Ql and M 2 > M 2 
11(b) : Ql > Ml > Ml 



tD 



tD 



r {mi,gi} 



Mi=0 



M,=0 



Mi=0 



r {mi,m 2 } 



Mi=0 



, T {mi,p 3 } 



Mi=0 

(3.39) 

with the shorthand notation 

J 3 D |m 1= = ^ ^ ^ °' ' ' M ' ' °) • (3 - 40) 

Again the hypergeometric functions collapse to Gauss' 2 Fi functions (see Eqs. ([A. 171) , 

( Pmp and qog )) 

F 2 (a, ^,^,7, 7', )2/ ) = 2 F 1 (a,/3 , ,7 , ,2/), 

5i (a, a', /?, 7, 6, 0, y) = 2 F X (a, a', 7, y) , 

H 2 (a, (3, 7, 5, e, 0, y) = 2 Fi (7, 5, 1 - a, -y) , 

and we recover, in region I, the result of Eq. (|3.16| ), and in region 11(b), the expected 
result (p^D for Mf < Q\. 

The M 2 -»■ limit: i/ 2 , i/ 3 ; Q?, 0, 0, M-j 2 , 0, 0) 

Taking the limit M 2 — > provides us with the integral relevant for the exchange of a 
heavy particle in the decay into two light particles. 

As usual, we could merely return to the system and, by setting m 2 = 0, solve it afresh: 
there are now (m + q + n) = 5 variables and still (n + 1) = 4 constraints leaving five 
single-sum solutions. However, it is simpler to take the M 2 - ► limit in the appropriate 
regions: 11(a) and 111(a), as can be seen from Eq. (|3.17Q . 

We then obtain 



11(a) 
111(a) 



Ql > Ml > M, 
Ml >Ql> M. 



I D 



tD 



M 2 =0 



M 2 =0 



j{mi,m 2 } _|_ j{m 2 ,p 2 } _|_ j{m 2 ,p 3 } 



M 2 =0 



r {m 2 ,pi} T {m 2 ,qi} 
3 ~r J 3 



M 2 =0 



where 



rD 



M 2 =0 



(3.41) 
(3.42) 



1° (p^p^-Ql 0,0, M?, 0,0). 

The hypergeometric functions collapse to 3 F 2 functions, according to Eqs. ( |A.10| ) 
and ( |Q3D 

Si (a, a', (3, 7, 5, x, 0) = 3 F 2 (a, a', (3, 7, 8, x) , 
S 2 (a, a', p, (3', 7, 0, y) = 3 F 2 (1 - 7, /3, /3', 1 - a, 1 - a', -y) , 
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and we obtain: 



if Q\ > Ml, region 11(a) 

4>i, Ql o, o, Ml o, o) = 4 mi} + 4 P2} + 4 P3} 



-i)* [Ql 



2 . ^-^-^-^ T f I/ X + Z/ 2 + I/ 3 - r - ^1 - V 2 ) V\^-Vx — Vz 



r (i/ 2 ) r (i/ 3 ) r (i? - i/i - 1/ 2 - ^3) 

/ DDL) Mf 

x 3 F 2 1 1/1, 1 + o- - D, a - —, 1 + uy + v 2 - —, 1 + v x + v 3 - — , 



1 J v y v y r fa) r(i*)r (§-1*) 

„ ( D D D Mf 

x 3^2 I ^2, — - ^3, 1 + ^2 - — , 1 + ^2 - ^3, 1 + y ~ v i ~ u 3, 

+ (-D* (-*! 



2\ 2 



r (r/ 2 - 1/3) r (f - v 2 ) 


r 






r> x )r> 2 )r 


(f 







X 3^2 f^'y _V 2, 1 + ^3 - y, 1 + ^3 + y -Vi - > ( 3 - 43 ) 

if Ml > Ql, region III (a) 

v 2 , v 3 ; Ql 0, 0, Ml 0, 0) = + 4 M 

( D D D Q 2 \ 

X 3 F 2 I Z/ 2 , Z/ 3 , vx + V 2 + V 3 ~ y > 1 + ^3 + ^2 - y , y , --r^ J 

r(, 2 + , 3 -f)r(f-, 2 )r(f-, 3 ) 



r (i/ 2 ) r (z/ 3 ) r (d - z/ 2 - z/ 3 ) 



( D D D Q\\ 

X 3^2 I y - v 2 , — - z/ 3 , D - v 2 - v 3 , 1 + y - i/ 2 - V 3 , -j^ I . (3.44) 

We can further check these results by setting one of the V{ — > to form bubble integrals 
or by taking one of the limits Mi — ► or Q\ — > 0. In each case, we recover the correct results 
presented in the earlier sections 



Discussion of system of constraints 

As we have already anticipated in Sec. [2.4. 2| , the homogeneous counterpart of the system of 
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constraints ( |3.3| ), with q 2 = q% = m 3 = 0, is not uniquely solvable. We can build the table of 
signs for qi and 



> 


< 


uncertain 


P2, qi, mi 


Pi, P3, m 2 




P2, P3, mi 


pi, qi 


m 2 


Pi, qi 


P2, P3, mi 


m 2 


Pi, P3, m 2 


P2, qi, mi 





but this does not determine the complete list of groups, even if some of them (first and 
second row, corresponding to region I and region 11(b)) are correctly predicted by the table 
(if we arbitrarily assume that m 2 in the second row has a positive sign). 



3.3 Equal-mass propagators: Mi = M 2 = M 3 = M 

In the special case of Mi = M 2 = M 3 = M we use the modified form given in Sec. |2.5|. For 



the most general case of unequal off-shell external legs, there are (n + q + 1) = 7 summation 
variables and (n+1) = 4 constraints. The template solution is easily obtained from Eq. fl2.58|) 
and the system of constraints is given by 

<?2 + <?3+Pl = -"I, 
<?l+<?3+P2 = ~V 2 , 

qi + q2+P3 = -V3, 
Pi + P2+P3 + qi + q2 + q3 = _ ? + m - ( 3 - 45 ) 



3.3.1 Ql = Ql = 

In this case, there are five summation variables (pi, p 2 , p%, qi and m), and four constraints, 
and the system admits four solutions out of the possible five. 

Based on two-particle cuts of the diagram, we see that there is a threshold at Q\ = 
4M 2 corresponding to producing propagators 2 and 3 on-shell. Of the four solutions, three 
converge when Q\ > 4M 2 , while the remaining solution converges when Ql < AM 2 . Recalling 
cr = vi + v 2 + z/ 3 (see Eq. ( |2.6| )), we have: 

if Ql > AM 2 

I°(u u u 2 , u 3 ; Ql, 0, 0, M 2 , M 2 , M 2 ) = / 3 {m} + + 
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*-r(<x-f)r( 



ui - u 2 ) r (f - v x 



T(u 2 )T(u 3 )T(D-a) 



( D 1 1 D D 4M r 

x 3^2 ( a - —, 1 + -{a - D), -(l + a -D),l + v 1 + u 2 -—,l + v 1 + v 3 - -, 



+ (-1)- 



r (1/3) r (z/i - 1/ 2 + ^3y 



(l. , 1, . D 

x 3 F 2 I z/ 2 , 1 + -(z/ 2 - ^1 - f 3 ), ^l 1 - ^1 + ^2 - fa), 1 + ^2 - ^3, 1 + y 



4M 2 
Of 



+ (-i)^(g? 



" 3 r-M 2 )T-i-^ 



r (i/ 2 - vz) r (^1 + v 2 - 
r (f 2 ) r iv x + V3) 



( 1 1 . D 

X 3 F 2 I I/3, 1 + -(1/3 - I/i - I/ 2 ), ^l 1 - ^1 + ^3 - ^2), 1 + Z/ 3 - Z/ 2 , 1 + y 



4M 2 



^1 - ^2, 



(3.46) 



if Q\ < AM 2 

I°(v h i, 2 ,v 3 ;Ql,a,0,M 2 ,M 2 ,M 2 



r'~ D 



_ An} 

( D a 1 + a Q 2 \ 

3 F 2 ( v 2 , 1/3, cr- — , -, — — , — — - I 



2 ' 2 



4M 2 / 



(3.47) 



where we have made use of formula Q2.64T ). This latter result agrees with that obtained by 
taking the limit Q^—^O, Q 2 — >0 in the general result given by Boos and Davydychev |13| 



for a triangle loop integral with three off-shell legs and a single mass M running round the 
loop. Equation ( |3.46| ) appears to be a new result. 



We note that, taking the limit v x — > in Eqs Q3.46|) and (|3.47|) , we reproduce the expected 
equal-mass bubble integral of Eqs. (|2.62 ) and (|2.63| ) 



^(0, V2, ^3; Ql 0, 0, M 2 , M 2 , M 2 ) = J 2 D (z/ 2 , 1/3; Ql M 2 , M 2 ), (3.48) 
while taking M — >0, only the first term in Eq. (|3.46 ) survives, yielding Eq. ([3 .11 ). 



However we observe that there are dangerous T functions in the second and third lines 
when z/2 = v 3 . Therefore to evaluate the integral when v 2 = v 3 = u, we introduce an 
additional regulator 5 such that v 2 = v + 5 and u 3 = v. As in the previous section, because 
the result does not depend on 5, the limit 5-^0 can be safely taken after the singularities 
have been canceled. 
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4 Conclusions 



Finally let us summarize what we have accomplished in this paper. Changing the number 
of dimensions D to evaluate loop integrals is well established and relies on the analytic 
properties of loop integrals. We have used this property to extend the number of dimensions 
to negative values as suggested by Halliday and Ricotta. As discussed at length in Sec. 



treating D as a negative even integer allows a multinomial expansion of the integrand in 
intermediate steps and, by expanding before and after loop integration we can identify the 
loop integral as an infinite series, together with constraints on the summation variables. The 
number of summation parameters is equal to the number of legs n plus the number of energy 
scales (m + q) in the loop, while there are n + 1 constraints, which can be read off from 
the Feynman graph. The form of the series is specified for arbitrary one-loop integrals and 
forms a template series into which specific solutions of the system of constraints are inserted. 
In this way, integration over the parameters is replaced with infinite sums. In each case 
we immediately identify generalised hypergeometric functions and show how to assemble 
the complete result valid in a particular kinematic region by considering the convergence 
properties of the hypergeometric functions. The procedure is as follows: 



1. Write down V, M. and Q of Eqs. (|2.9|) , (12.101) and ( |2.11|) by inspection of the graph. 



2. Write down the template solution for the particular process. In other words, copy out 
Eq. Q2.23 ) inserting the correct mass scales and number of terms in V, A4 and Q. 



3. Construct the system of constraints by counting powers of Xi in Eqs. ( |2.20| ) and ( |2.21 



4. Solve the system of constraints and insert each solution into the template solution, one 
at a time. 

5. Construct Pochhammer symbols and identify the generalised hypergeometric function. 

6. Flip all of the T functions in prefactor according to Eq. (|2.33|). 



7. Group the solutions according to their regions of convergence. 

8. Evaluate the hypergeometric functions for the specific parameters of interest. 



Steps 1-6 are very straightforward and easily achieved with a computer program. Step 7 
requires a little more thought, though for the cases we have studied here the convergence 
regions were easy to identify. To make the procedure useful for phenomenological studies, 
it is necessary to evaluate the hypergeometric functions for specific values of the parame- 
ters. In particular, an integral representation is required. The mathematical literature for 
Euler integral representations of hypergeometric functions with more than 2 or 3 variables 
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is quite sparse, and it may be necessary to select a complex integral representation for more 
complicated functions. 



More interesting is the application of NDIM to integrals with more than one loop. Suzuki 



and Schmidt |23|, |24j] have made some steps in this direction, although the integrals they have 
considered are largely of the one-loop insertion type. Given that quite powerful results for 
one-loop integrals are achieved so easily and generally, we expect that NDIM can play a 
role in simplifying the task of calculating two- (or more) loop integrals that are necessary to 
make more precise perturbative predictions within the Standard Model. 



Acknowledgements 

We thank J. V. Armitage for stimulating discussions and insights into the field of hyperge- 
ometric functions and P. Watson and M. Zimmer for useful conversations. We acknowledge 
the assistance of D. Broadhurst and A. Davydychev regarding generalised hypergeomet- 
ric functions. C.A. acknowledges the financial support of the Greek government and CO. 
acknowledges the financial support of the INFN. 



A Hypergeometric definitions and identities 

The purpose of this appendix is to give sufficient information to evaluate the general loop 
integrals presented in Sees. |27| and 0. In Sec. |A.1| we give the definitions of the hypergeomet- 



ric functions as a series together with their regions of convergence. Integral representations 
are provided in Sec. |A.2j together with a description of how to evaluate the integrals in the 
general case. For specific choices of the z/j the general hypergeometric functions often simplify 
and some useful identities and analytic-continuation formulae are collected in Sec. [A.3| . 



A.l Series representations 

The hypergeometric functions of one variable are sums of Pochhammer symbols over a single 
summation parameter m, like, for example, 

n / a \ (a,m)((3,m) x m ,. 

2Fl(Q ' A7 ' l) = g„ ( 7 'm) {AA) 

m =o W,m) (7',m) ml 

which are convergent when |x| < 1. 
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We also meet hypergeometric functions of two variables which can be written as sums 
over the integers to and n: Fi, % = 1, . . . , 4 are the Appell functions, H 2 a Horn function and 
Si and S 2 generalised Kampe de Feriet functions: 

(a, to + n)((3, m)((3', n) x m y n 

to! n! 



F 1 {a,p,(3',y,x,y 

F 3 (a,a',f3,f3',>y,x,y 
F 4 (a,P,y,Y,x,y 
H 2 (a,/3,y,y',8,x,y 
St {a,ot,(3,y,8,x,y 

S 2 (a,a',P,P',y,x,y 



n (7, m + n 

m,n=0 \ I ' 

y, (a, to + n)((3,m)((3' ,n) x m y 

m,n=0 



(7, to) (7', n) to! n! 

'a,m)(a',n)(j3,m)((3',n) x m y n 
(7, m + n) to! n! 

a, to + n) (/?, m + n) x m y n 



E 

m,n=0 
00 

y ■ 

m % (7,m)(Y,n) to! n! 

x (a, m — n)(/3, to) (7, n)(7', n) x m y n 

to! n! 



E ■ 

m,n=0 



E 

m,n=0 
oo 

E ■ 

m,n=0 



(5, to) 

Vy, TO + n)(o; / , m + n)(/3, to) x m y 71 
(7, to + n) (5, to) to! n\ 

a,m — n) (a', m — n) {(3, n) (/?', n) x m y n 



to! n! 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



(7, to — n) 

These series converge according to the criteria collected in Table [l], The domain of con- 



Function 


Convergence criteria 


F x , F 3 


\x\ < 1, \y\ < 1 


F 2 , S 1 


\x\ + \y\ < 1 


F 4 


yN + y\y~\ < 1 


H 2 , S 2 


-\x\ + l/\y\ > 1, \x\ < 1, |y| < 1 



Table 1: Convergence regions for some hypergeometric functions of two variables, 
vergence of the Appell and Horn functions are well known. That one for Si and S 2 may be 



worked out using Horns general theory of convergence [28 



When one of the arguments vanishes, then the hypergeometric function collapses in a 
straightforward way. For example, if y = in Eq. ( |A.8[) , then only the first term of the series 
in n contributes and we are left with the relation, 



Si (a, a', (3, 7, 5, x, 0) = 3 F 2 (a, a', (3, 7, 5, x) 

Similarly, we have 

Si (a, a',/?, 7, 5,0, y) = 2 F X (a, a', 7, y) 



(A.10) 
(A.11) 
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S 2 (a,a',P,y,^,x,0) 


= 2-^1 


[a, a', 7, x) 


(A.12) 


S 2 (a,a',p,p',y,0,y) 


= 3F2 


(1-7,/3,/3',1-a, l-a',-y) 


(A.13) 


(a, A ^,7,^0) 


= 2-F 


[<x,P, n/,x) 


(A. 14) 


F\ (a, /3, /?', 7, 0, y) 


= 2-^1 


(a,/?', 7,1/) 


(A.15) 




= 2-^1 


(a,/3,7,x) 


(A.16) 




= 2-^1 


(a,P',j',y) 


(A.17) 


F 3 (a,a',P,P',rf,x,0) 


= 2F1 


(a,P,^,x) 


(A.18) 


F 3 (a, a', p, (3', 7, 0, y) 


= 2F1 


(u',P',i,y) 


(A.19) 


F 4 (a, /?, 7, 7', 0) 


= 2F1 


(a,/?,7,x) 


(A.20) 


F 4 (a, /3, 7, 7', 0, y) 


= 2F1 


(a, 7', 2/) 


(A.21) 


#2 (a,/3,7,7',5,x,0) 


= 2F1 


[a, P, 5, x) 


(A.22) 


#2 (a, P, 7, 7', 5, 0, y) 


= 2F1 


(7,7', 1 - a, -y) . 


(A.23) 



When one of the parameters vanishes producing a Pochhammer (0, n), then the series in 
n also terminates. If we have (0,m + n) then both series terminate. For example 

F 2 (a,0,P', 7 ,y,x,y) = 2 F, (a, P' ,7' ,y) (A.24) 
F 1 (0,P,P', 7 ,x,y) = 1. (A.25) 



A. 2 Integral representations 

Euler integral representations of the hypergeometric series of one and two variables are well 
known and we list them here. We know of no integral representation for the H 2 

function and for the closely related S 2 function. 

2 F X {a, p, 7, x) = — ^ x C duu^-\l - uy- p -\l - uxY a 



3 F 2 (a,P,y,6,e,x 



T(P)T( 1 -P) Jo 

Re(p) > 0, Re(7 - p) > 0. (A.26) 

r (g r (e) 
r(/3)r(5-/3)r( 7 )r( e -7) 

x / du f dvu p - 1 v^-\l-v) s - p - x [l-v) e -^- x [l-uvx)- a 
Jo Jo 

Re(p) > 0, Re(S - p) > 0, Reft) > 0, Re(e - 7) > 0. (A.27) 
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F 1 (a, (3, p, 7, x, y) = — ^ /' du u"-\\ - u)t-°-\\ - ux)^{\ - uy)^' 

1 (a 1 (7 — a) Jo 



F 2 (a, 0, ft, 7, y, x, y) 



r (a) r (7 - a) 

Re(a) > 0, Re(7 - a) > 0. (A.28) 

E (7) r (y ) 
r (p) r 09') r ( 7 - /?) r (y - 



X 



F 3 (a,a',P,P',-f,x,y) 



[ du [ dvuP-W-^l-uy-P-^l-vy'-P'-^l-ux-vy)-" 
Jo Jo 

Re((3) > 0, Re((3') > 0, Re( 7 - (3) > 0, Re( 7 ' - > 0. (A.29) 

r(7) 



r (0) T (/?') T (7 - - /?') 

fU+V<l 



X 



// " duduu^-V-^l-u-u) 7 -^- 1 ^-^)-^!-^)- 

J Ju>0,v>0 

Re(/3) > 0, Re(/3') > 0, Re( 7 — /? — /?') > 0. (A.30) 

F 4 (a, A 7, 7', *(1 - y), y(l - x)) = r ( a ) r (^7- a)V (Y - /?) 

x /"^u f 1 dvu^v 13 ' 1 ^ -u)^ a -\l -v) 1 '-^ 1 
Jo Jo 

X (1 - - Vyf-^' +1 {l — MI — vy )7+7'-a-^'-l 

Re(a) > 0, Re(/3) > 0, Re(7 - a) > 0, Re(Y - /3) > 0. (A.31) 

S 1 (a,a',P, 1 ,5,x,y)= \" / duu a ~\l - U y- a - 1 F 2 (a' : 0, 1, 5, 1, ux, uy) 

1 (a J 1 (7 — a ) Jo 

r ( 7 ) r (5) 



r»r( 7 -a)r(/3)r(5-/3) 

fl /•! 



x [ du [ dvu a - 1 vP- 1 (l-uy- a - 1 (l-v) s - l3 - 1 (l-uvx-uy)- a ' 
Jo Jo 

Re(a) > 0, Re(7 - a) > 0, Re(/3) > 0, Re(5 - /3) > 0. (A.32) 
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A. 2.1 Example of explicit evaluation of an integral representation 



In working out the integral representation for hypergeometric functions in D = 4 — 2e 
dimensions, we have often to deal with the e expansion of integrals of the form 



I(x) 
d{u) 



u 



du d{u) f(u), 
1+Q£ (1 - u)- l+Pt 



(A.33) 
(A.34) 



where a and j3 are real numbers and f(u) is a smooth function in the domain < u < 1: in 
particular, it is finite at the boundary points. 

The procedure to deal with this kind of integrals is quite standard. The integral has a 
pole in e when the integration variable u approaches either of the end points. We concentrate 
first on the point u = 0, and we rewrite the integral in such a way to expose the pole in e 



(A.35) 



I(X) = j dud(u) f(0) + J Q dud{u) [f(u) - /(0)j = /[!] + I m . 

The integral J[i] can be easily done 

T(ae)T((3e) f(0) a + (3 T (1 + a e) T (1 + (3 e) 



f(oy- 



(A.36) 



r((a + P)e) e a(3 T (1 + (a + (5) e) ' 

and the integrand of i™ is now finite in the limit u^O. In fact, we can make a Taylor 
expansion 

f(u) - /(0) = u/'(0) + |T(0) + . . . = ug(u), (A.37) 

and write ip] as 



dud(u) ug(u) 



duu ae (l-u) 



9W 



(A.38) 



We repeat now the same steps done for Eq. ( |A.35|) with respect to the point u — 1, to obtain 



I [2] = I duu a \l-u)- 1+ ^g{\)+ j duu a \l-uY 1+ ^ g{u)-g{\) =I n +I W (A.39) 



The integral J[3] gives 



= r(l + ae)r(/36) /(l)-/(0) r(l + a6)r(l + /3e) 
[3] 9{ ' r(l + (a + /3)e) /3e r(l + (a + /3)e) 



(A.40) 



while J[4] is finite at u — >• 1 



dwW ae (l-w)^ e /i(M), </(«) = (l-u)/l(w), 



(A.41) 
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and can be solved with an e expansion of the integrand. Adding all the contributions together 
we have 



I(x) 



where 



1 ^ ( o) +Q /(i)] r(1+a£)r(1+/ * £) ■ r l 



a (3 e 



V (l + (a + P) e) 



+ / duu ae {l-uf e h{u), (A.42) 
Jo 



h(u) 



) = 7T~ r (/(«) - (1 - u) f(0) - u /(l) 
u(l — u) \ 



(A.43) 



In the case where we have two integration variables, the procedure outlined above can be 
re-iterated in a straightforward manner. To illustrate the procedure, we evaluate explicitly 
the F 2 functions of Eq. ( ggg ) to O (e 2 ). 

The integral representation for F 2 (see Eq. ( |A.29 )) is given by 

e 2 r(l- e ) 



F 2 (1, l,e,e + 1, 1 - e,x,y) 



r (1 + e) T (1 - 2e) 



where 



(A.44) 
(A.45) 



and 



I(x,y) = / dudv d(u,v) f(u,v), 
Jo 

d(u,v) = v- 1+e (l-u)- 1+e (l-v)- 2e 
f(u,v) = (l-ux-vy)-\ 

and I(x,y) must be computed to 0(e°). In order to expose the poles (see Eq. (|A.35|) ), we 
add and subtract the value of the finite function f(u,v), computed at the boundary points, 
in the following way: 

I(x,y) = f 1 dudvd(u,v){[f(l,0)] + [f{u,0)-f(l,0)] + [f{l,v)-f(l,0) 



+ 



/(«,«)- /(«, o)-/(M) + /(i,o)]} 

= / [1 ]+/ [2 ]+/[3] + /[4]. (A.46) 

We are now in a position to evaluate the single contributions in the square brackets. In fact 



'[i] 



42] 



'[3] 



4i] 



(1 -xy 1 / du{l-u)~ 



■l+e 



dvv- 1+e (1 - v)- 2e = (1-x) 



. 1 r(l + e)r(l-2e) 
e 2 T(l -e) 



-x r (1 + e) T (1 - 2e) /• 1 ^(l-w) f 



1 — x eT(l — e) jo (1 — ux) 

(1-a;)- 1 /-i , v e (l-v)- 2 " 
dv 



xy 
1 — x Jo 



1 — x — v y 



du dv (1 — u) e v e (1 — v) 



-2e 



(vy + ux + x — 2) 



(1 — ux)(l — x — vy)(l — v y — ux) 
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(A.47) 



The remaining integrals are finite in the limit e — > 0, so that we can make a Taylor expansion 
to O (e) for the integrands of and i"^ , and we can put directly e = in Iw . Recalling the 
definition of the dilogarithm function 

Li 2 (x) = - f dz lQg(1 - Z) x<l, (A.48) 
Jo z 

it is straightforward to carry on the last integrations and express the result in terms of Li 2 
functions, as done in Eq. (|3.8 1|). 



A. 3 Identities amongst the hypergeometric functions 

There are three kinds of identities that relate hypergeometric functions. First there are 
analytic continuations which connect functions in different regions of convergence. Second 
are reduction formula which allow the functions to be expressed as simpler series for certain 
values of the parameters. Finally there are transformations which relate the same functions 
with different arguments. 



A. 3.1 Analytic continuation formulae 

Here we give only those analytic continuation properties that relate the argument and inverse 
argument. Gauss' hypergeometric function has the following analytic continuation properties 
(see for example [p0| ) 



2 Fx (a, (3, 7 , z) = {-z)- a V J?)?Jf " a*i fa, 1 + a - 7, 1 + a - (3, - 

r (p) r (7 - a) v z 

\axg(-z)\ < Ti, (A.49) 

2 F X (a, (3, 7 , z) = z - ^(7)r(7-a-/3) p f 1 + a _ 1 + a + p _ 

r (7 — a) r (7 — p) \ z 

. a _ 7 n , 7 -a-p r (7) r + g - 7) / 1 

+2 7 1-2 w > wm 2-Fi 7-a, 1+y-a-P, 1-- 

r (a) r (/?) v 2: 

|arg(^)| < 7r, |arg(l -2)] < 7r. (A.50) 

The corresponding analytic continuation of the hypergeometric functions with two vari- 
ables are summarised in Table [| There are many possible analytic continuations; however, 
we list only those that are relevant to link the groups of solutions for the two- mass triangle 
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integral discussed in Sec. |37^, that is the connections between the Appell and Horn functions. 
The others are easily derived by summing the series with respect to one of the summation 
variables to obtain an 2 F^, applying Eq. ( A. 49 ), rewriting Gauss' hypergeometric function 
as a series and reidentifying the double series. We see that these functions appear to form a 
group. 



Function 


Continued in terms of 


F 4 (x,y) 


F 4 (x/y, 1/y), F 4 (y/x, l/x) 


F 3 (x,y) 


H 2 (l/x, -y), H 2 (l/y, -x), F 2 (l/x, 1/y) 


H 2 (x,y) 


F 2 (x, -1/y) 


F 2 (x,y) 


Si(-y/x, l/x), H 2 (y, -l/x), St(-x/y, 1/y), H 2 (x, -1/y) 


H 2 (x,y) 


F 3 (l/x, -y), S 2 (l/x, -xy) 


Si(x,y) 


F 2 (-x/y, 1/y) 




H 2 (l/x, -xy) 



Table 2: Analytic continuation for the hypergeometric functions of two variables. 



F 4 (a,(3,<y,<y',x,y) 



r ( 7 ') r {(3 - a) 



+ 



F 3 (a,a',(3,{3',y,x,y) 



+ 



r(y 


-a)T{(3) 


r ( 7 ') T(a-P) 


r(y 


-P)T(a) 


r(/3 


-a)r( 7 ) 


r( 7 - 


-a)T((3) [ 




-/?)r( 7 ) ( 



(-y)~ a F 4 [a, a + 1 - j', j, a + 1 - (3, 



x 1 

y y y 

X l" 



-y)- p Ft /?, (3 + 1 - i, 7, P + 1 - a, - - (A.51) 

V y y 



—x) a H 2 ( a + 1 — 7, a, a', (3', a + 1 — (3, — , —y 



~ -3)T(a) { - x) ~" H2 { P + 1 " 7 ' A a '> 13 



+ 1 - a 



x 

,-,-y 

X 



(A.52) 



H 2 (a, (3, 7 , i, 5, x, y) = ^ j) T (1 a > (y)^F 2 ( a + % (3, j, 5, 7 + 1 - Y, x, --) 

r(i-«-7)r(7') v y j 

+ In' 1 '^^^ '"'^ ( a + V, & 7', S, i + 1 - 7, x, --) 

r(i-a-7')r(7) v yJ 

(A.53) 

F 2 {a, (3, P>, 7, 7', x, y) = * ^ ~ (-*)-*& (a, a + 1 - 7, /?', a + 1 - (3, 7', -|, I 

+ rl^-^rS ^"^ 2 ( a - * * * ^ + 1 - 7 ' V ' y ' 4) (A - 54) 
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H 2 (a, (3, 7, y, 5, x, y) = ^ — ^tftIt (~ x ) f «, a + 1 - 5, 7', 7, a + 1 - /3, -, -xy 

r d— a r « V x 



r(* 


-a)r(/3) ( 


r( a 


-)9)r(5) 


r(* 


-/?)r» l 


1 (a 


- a) r (7) 


r( 7 


- a) r (a') 




- «') r (7) 


r( 7 


- a') T (a) 


r(a' 


-a)r( 7 ) 


r( 7 


- a) r (a') 


r(a 


- «') r ( 7 ) 



(A.55) 



Sx(a,a',(3,"{,5,x,y) = ^ / ' { (-y) a F 2 ( a,/?, a + 1 - 7, 5, a + 1 - a', --, - ] 

1 V 2/2// 

+ ;r"^!'l N^ f ft a' + 1 - 7, 5, a' + 1 - a, --, - 

I y y 



(A.56) 



S 2 (a, a', (3, (3', 7, x, y) = -rwr a (~ x ) "-^ 2 ( a, a + 1 — 7, /?', /3, a + 1 — 0/, -, -xy 



T (7 — a') r (a) V x J 

(A.57) 



A. 3. 2 Reduction formulae 

The F4 functions describing the massive bubble and the off-shell massless triangle have the 



following reduction formulae which leave a single remaining Euler integral at most |3(| |3T 

F 4 a, (3,^, (3 





X 






2/ 




(1- 


x)(l 


-vY 


(i- 


-x)(l 


-2/) 






= (1 


— x] 


«(!_ 


2/) a * 




X 






1/ 




(1- 


x)(l 

X 


-vY 


(1- 


-x)(l 
2/ 


-2/) 


(1- 


x)(l 


-vY 


(1- 


-x)(l 


-2/) 



^4 (a, ft a, ft — - , — J„ - ) = (1 - xy)- x (l - a^(l - y) a , (A.59) 

F 4 (a,ftftft 



F 4 a,ftl + a-ftft 



(1 - x) Q (l - y) a 2 F l (a,l + a-(3, (3, xy) , (A.60) 
x y \ 



(l-x)(l-yY (l-x)(l-y)J 

= (1 - 2/) a 2 Fi U(3,l + a- (3, -^^) ■ (AM) 
Similar reductions for the other functions of two variables are 

F x (a, (3, (3', (3 + (3', x, y) = (1 - y)-^ (a, (3,(3 + (3', (A.62) 
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F 2 (a, ft ft, 7, a, x, y) = (1 - y)^' F 1 (ft a - ft, ft, 7, x, (A.63) 

F 2 (a, ft ft, a, 7', s, y) = (1 - x)^F 1 (ft, (3, a - (3, 7', j^,y) (A. 64) 

F 2 (a, (3, ft, (3, 7', x, y) = (l- x)- a 2 F l (a, ft, i, ^) (A.65) 

F 2 (a, ft ft, a, a, x, y) = (1 - x)^(l - y)^' 2 F 1 fft ft, a, ^—^L-—^j (A.66) 

F 2 (a, ft ft, a, ft, x, y) = (1 - - x - y)"' 3 (A.67) 

F 2 (a, ft ft, ft ft, x,y) = (l-x- y)~ a (A.68) 

F 3 (a, 7 - a, ft ft, 7, x, y) = (1 - y^Fi (a, ft ft, 7, x, (A.69) 

F 3 (a, 7 - a, ft 7 - ft 7, x, y) = (1 - yr+^Fx (a, ft 7, x + y - xy) , (A.70) 
while for certain values of the parameters the H 2 function reduces to an F 2 or F\ 

H 2 {a, /3, 7, 5 - a, 5, x, y) = (1 - x)~ p F 2 (s - a, ft 7, 5, 1 - a, , -y) , (A.71) 

H 2 (a,f3,7,l-a,5,x,y) = (1 + y)^F 1 (ft a,<y,5,x, . (A.72) 

The S*i and 5*2 functions we have introduced are less well known. From the integral 
representation or manipulating the series using Eq. ( |A.76j ) we find 



Si (a,a',@,a,5,x,y) = (1 - y) a ' 2 i<i \a r , (3, 5, > (A.73) 

St (a, a', ft 7, ft x, y) = 2 Fi (a, a', 7, x + y) , (A.74) 
S 2 (a, a', ft ft, a,x,y) = (1 - x)~ a ' 2 F 1 (ft ft, 1 - a', -y(l - x)) . (A.75) 

A. 3. 3 Transformation formulae 

A useful formula connecting Gauss' hypergeometric function to itself is 

2F1 (a, /3, 7, z) = (1 - zY a 2 F X fa, 7-/3, 7, 

= (1-2)-" 2 Fi( T -a, ft 7, : 



2-1 

(l-zr^ 2 F 1 (7-«,7-A7^). (A.76) 
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Using the above results, then 

^1 (a,(3,(3',l,x,y) = (l-^)- /3 (l-l/)- /3 'F 1 f7-«,/9,/9',7,^ T ,^ T N ) 

\ x - 1 y - 1J 

\ x — 1 x — \, 



= (1 - y)~ a F 1 [ a,/3,7 - /3 - /3',7, - — t, — ] (A.77) 



F 2 (a,P,P' nn ',x,y) = (l- x )-F 2 (a, 1 -P,P' nn ',^ 1 , T ^ 

= (l-y)- a F 2 (a,p,i-p' nn ',-^,-^- 
\ \-y y-\ 



[l-x-y)-°F 2 [a, 7 -/?,-/ -^,7,-/, ' " 



x + y — 1 x + y — 1 J 
(A.78) 
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